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Abstract 



We relate Bondi systems near space-like infinity to another type of gauge conditions. While 
the former are based on null infinity, the latter are defined in terms of Einstein propagation, 
the conformal structure, and data on some Cauchy hypersurface. For a certain class of time 
symmetric space-times we study an expansion which allows us to determine the behavior of 
various fields arising in Bondi systems in the region of space-time where null infinity touches 

^ ■ space-like infinity. The coefficients of these expansions can be read off from the initial data. 

We obtain in particular expressions for the constants discovered by Newman and Penrose 

5h . (NP-constants) in terms of the initial data. For this purpose we calculate a certain expansion 

^W' introduced in [k|| up to 3rd order. 



1 INTRODUCTION 2 

1 Introduction 

Most studies of gravitational fields near null infinity are based on the use of "Bondi-type" coordi- 
nates. In the first investigations of the behavior of the field near null infinity (cf. |[|], ||ll|, P4 ) 
Bondi-type coordinates played a crucial role in the specification of the fall-off behavior of the field. 
The characterization of the asymptotic behavior of gravitational fields near null infinity in terms of 
the conformal geometry subsequently suggested by Penrose (iQ, Ijl^l) does not require the use of 
such a specific class of coordinates. Nevertheless, Bondi-type coordinates are usually also employed 
in this context because they allow us to exploit in a convenient way certain features of the null cone 
structure. If the gravitational field is, however, to be analyzed in detail in the region where future 
and past null infinity J^^ "touch" space-like infinity, and if this is to be done such that J'~ and J'^ 
are treated on an equal footing, Bondi-type coordinates are not particularly helpful. Already in 
the simplest non-trivial case, that of the Schwarzschild solution, the use of double null coordinates 
leads to difficulties. 

In iQl an initial value problem for the conformal vacuum field equations has been formulated 
which is designed to analyze near space- like and null infinity the Einstein propagation of asymptot- 
ically fiat data on a Cauchy hypersurface 5 in a finite picture. In this setting, which is based on cer- 
tain conformally invariant structures, space-Hke infinity is represented by a cylinder J ~] — 1, l[xS'^ 
such that the sets J^ ~ K x S^, representing future resp. past null infinity, "touch" the cylinder 
at its two boundary component J^ = {il} ^ S^- Though the underlying facts about the evolution 
equations which have been used here hold for much more general situations, the picture has been 
analyzed so far under certain simplifying assumptions on the initial data. The data are assumed to 
be time-symmetric and the conformal structure, which then represents the free datum, is assumed 
to extend smoothly through space-Hke infinity such that the latter is represented by a point i in an 
extended manifold 5" = 5 U {i}. The cyHnder / is obtained by blowing up the point i to a sphere 
J° ~ {0} X S^ and by smoothly extending the solution in a particular geometric gauge. 

It can be seen already under these assumptions on the data that the new picture allows us to 
relate near /=*= properties of the data on S, which touches J at /", to properties of the field on null 
infinity by solving a hierarchy of differential equations on /. These equations have been used in |Q| 
to derive certain "asymptotic regularity conditions" for the initial data whose imposition prevents 
a certain class of logarithmic singularities of the field at the sets /* from arising. However, it 
still has to be shown that the asymptotic regularity conditions ensure a time evolution of the data 
which extends near space-like infinity smoothly to null infinity. 

In the present article we analyze the consistency of the early investigations of fields near null 
infinity with the picture developed in [|j and we demonstrate to some extent the efficiency of the 
latter in calculating near space-like infinity quantities on null infinity from the given data. For this 
purpose we make two different types of assumptions. On the one hand we shall consider space- 
times arising from time symmetric vacuum data as described above which satisfy the asymptotic 
regularity conditions. Our calculations of fields on the cylinder / rely only on these assumptions. 
On the other hand we shall assume that these data develop into solutions which admit a smooth 
conformal structure at null infinity and that the gauge conditions proposed in [Q extend in a 
smooth and regular way to J^^ . We expect that our analysis will contribute information on the 
solution process which in the end will allow us to remove the second type of assumptions and to 
show that the existence of the smooth evolution can be derived solely from assumptions on the 
initial data. 

The present article can be divided into three different, though related, parts. 
— In yi an expansion of the field near space-like infinity in terms of a "radial" coordinate p, which 
vanishes on the cylinder / representing space-like infinity, has been introduced. We calculate the 
coefficients of this expansion to third order. This calculation is not only of interest because it 
allows us to study the NP-constants, which will be discussed below, but also because it provides 
some information on the smoothness of the evolution near null infinity for fields arising from 
data subject only to our first type of assumptions. Though the asymptotic regularity conditions 
referred to above exclude certain types of logarithmic singularities in the evolution near /, there 
exists another potential source of singularities. To show that in fact no further singularities can 
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arise at any order, it is clearly of interest to understand the situation for the first few orders of 
the expansion. The potential singularities should show up for the first time at the order of our 
calculation. Our calculations show that at this order they are in fact excluded by the asymptotic 
regularity conditions. 

We note that our expansion of the field near space-like and null infinity, which we carry out 
in terms of the conformally rescaled fields and associated gauge conditions, can be translated 
into an expansion of the field near space-like infinity in terms of the "physical" field and suitable 
coordinates. We shall not carry out such a translation because the main point of our consideration 
is the fact that we can relate quantities on null infinity to the data on S. 

— Bondi-type coordinates and certain related frame fields (cf. the definition of the "NP-gauge" 
below) are based on the structure of null infinity. The gauge conditions in [Q (cf. the definition 
of the "F-gauge" below) are based on Cauchy data, the Einstein equations, and certain properties 
of conformal structures. We discuss in general terms how to construct near null infinity the 
transformation from the F-gauge into the NP-gauge. Using the expansion referred to above we then 
obtain expansions near /"^ of various quantities given in the NP-gauge in terms of the coordinates 
arising in the F-gauge and coefficients which are given directly in terms of the initial data on S. 
We note that these expansions imply expansions of quantities of physical interest on null infinity 
such as the Bondi-energy- momentum, the angular momentum (cf. |[l9J for various suggestions), 
the radiation field, etc. in terms of the coordinate p on null infinity, which vanishes at J"^, and 
coefficients derived from the initial data. 

Since we need for our considerations quite detailed information on the structure of the initial 
data near space-like infinity, our explicit calculation are done only for time-symmetric data. How- 
ever, many of our considerations apply also to more general situations and as soon as sufficient 
information on data with non- vanishing extrinsic curvature becomes available (cf. [^), we shall be 
able to derive by similar calculations relations between fields on J^~ and J^ . These relations will 
contain non-trivial information on the evolution process. 

— As a specific appHcation of this discussion we reconsider the constants which have been associated 
by Newman and Penrose with asymptotically simple space-times (cf. ||12|, [Q). The NP-constants 
are given by certain integrals over spherical cuts of null infinity and have been shown to be ab- 
solutely conserved in the sense of being independent of the choice of cut. We derive for them 
expressions in terms of the initial data on S. Such expressions have been given already in the 
static case in [Q. We derive analogous expressions for a much more general class of space-times 
arising from time-symmetric initial data. For these data the time evolution of the field is in general 
not known explicitly as it is the case in the presence of a time-Hke Killing vector field. The fact 
that we can nevertheless obtain expressions in terms of the data illustrates to some extent the 
efficiency of the new picture. Though various authors (cf. |g|, ^^, ||2J) discuss these constants 
from different points of view, no consensus has been found concerning their geometrical/physical 
significance. Whether our discussion will help clarify the meaning of the NP-constants remains to 
be seen. One of our main reasons for looking at them is the expectation that they may play a role 
in the construction of space-times. In numerical calculations they may certainly provide a check 
on the numerical accuracy. 

2 Relating different gauge conditions near null infinity 

We begin by giving an outline of the finite, regular initial value problem near space-like infinity. 
This has been introduced in the article y, to which we refer for more details. It involves a gauge 
which we refer to as the F-gauge. We then recall the NP-gauge, employed in jljl, to discuss the 
gravitational field near null infinity. Finally, we discuss how the NP-gauge is related to the F-gauge. 

2.1 The regular finite initial value problem near space-like infinity 

We want to discuss asymptotically fiat solutions {M,g) to Einstein's field equations i?^,y = 
in a neighborhood Ma of space-like infinity which covers parts of future and past null infinity. 
The solutions arise from asymptotically fiat data on a smooth space-Hke Cauchy hypersurface 
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S C M which are such that the intrinsic conformal structure on S admits an extension with 
a certain smoothness to a smooth compact manifold S obtained from S by adjoining a point i 
which represents space-hke infinity, S = S U {i}. We assume that the solution, i.e. the evolution 
in time of these data, possesses a smooth conformal extension (M, g, Q) such that we can write 
M ~ M U J^ U >/+, where J^ ~ M x 5^ represent future respectively past null infinity and 9 
denotes a smooth "conformal factor" on M such that O > and g — Q^ g on M while = 0, 
de^Oon J^±. 

To analyze in detail the consequences of the field equations in a neighborhood of space-like 
infinity which covers parts of J^^ , the situation above has been discussed in the |Q| in terms of a 
certain principal fiber bundle M^ -^ Ma with projection tt, 4— dimensional base space Ala, and 
bundle space M^ which is a 5— dimensional manifold with boundary and edges. To describe this 
setting further we need to introduce some notation. 

We employ the two-components spinor and space-spinor formalisms as used in Q where Cah, 
e°^ are the antisymmetric spinors with eoi = 1, e°^ = 1. We set r°° = eo" eo' ° + ei ° ei' " . By 
SU{2) will be denoted the group of 2 x 2 matrices t — (t"- f,) satisfying 

^acf" bf^ d^ ^bd, Tact"" bf^ d = Tbd, 

and by C/(l) its subgroup of diagonal matrices. A basis of the Lie-algebra of SU{2) is then given 
by the 2x2 matrices 

"1 = K ^ ) ' "' ^ K 1 0^ ) ' "' ^ K -* ) ' ^^-^^ 

of which Us generates U{1). 

In the following will be described in detail the regular finite initial value problem at space-like 
infinity formulated in |Q . Though we shall remark in passing on the construction of the manifold 
Ma and the underlying gauge conditions, we refer for the full details to the original article. The 
manifold M^ is given by 

M'a = {{T,p,t) e M X R X SU(2)\ < p < a, -- < T < -}, 

P P 

where a is a positive real number and uj = ^{p,i) a smooth non-negative function, given below, 
such that ^ extends to a smooth positive function with ^ ^ 1 as p ^ 0. By p and r will also 
be denoted the projections of M'^ onto the first respectively second component of IR x M x SU{2). 
Then any coordinate system on SU{2) will define together with the functions p and r a coordinate 
system on M^. There will, however, arise no need for us to introduce coordinates on SU{2). We 
denote the projection onto the third component of M x M x SU{2) by t and regard the S'[/(2)-valued 
function i as a "coordinate" on M^. 

The natural action on the right of U{\) on SU{2) induces a smooth action of U{1) on Af^. The 
quotient M^/[/(l) under this action will be denoted by Ma and the induced projection of M^ onto 
Ma by TT. We shall write N = 7r(iV') for any subset TV' of M^. The following subsets of M^ will 
be important for us: 

j'^ = {t = ±-,p>0}~Mx S'^ 

p 

I' = {\t\ < 1, p = 0} ~ M X S•^ /'± = {r = ± 1, p = 0} ~ S•^ 

C' ^{t = 0}, /° = {r = 0, p = 0} = C" n /' ~ 5*^ 

Because they cover only a part of null infinity close to space-like infinity, we should have denoted 
the first sets more precisely by J7^* but we dropped the subscript a for convenience. By definition 
the part of the physical manifold M which is covered by Ma is given by Ma = Ma \ [J^ U J^ U / U 
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/^ U /+), the sets J^^ represent future resp. past null infinity while the set / represents space-like 
infinity for Ma and the metric induced on it by g. Thus Ma covers a neighborhood of space-Uke 
and null infinity in M. The edges /* ~ 5*^ of Ala at which future resp. past null infinity touches 
space-like infinity will play an important role in the following. We shall refer to the set C as the 
initial hypersurface since by definition C n Ma = C \ l'^ = S O Ma. There exists a neighborhood 
Ba of i in S' and smooth surjective map n' : C ^ Ba which is injective on C \ /'^ and which maps 
/° onto i. 

As described in |Q|, the manifold Af^ is obtained essentially by lifting Ma into the bundle of 
normalized (with respect to Cab) spin frames. The set / ° ~ SU{2) corresponds to the set of 
normalized (with respect to Cat and Tab) spin frames at the point i. With each such spin frame we 
associate a unit tangent vector of 5 at i. With this vector we associate in turn a curve through i 
in Ba and extend the spin frame along this curve by a certain transport process. Thus we obtain 
spin frames at each point of Ba \ {i}. These frames are transported off Ba \ {i} ~ C \ /" into the 
space-time Ma by a certain propagation law along conformal geodesies orthogonal to C. The latter 
are given in our description of M^ by the curves p = const., t = const, with t a natural parameter 
along them. Since for given unit tangent vector at i the spin frame defining it is determined up to 
a phase factor, the spin frames at points of Ma \ (/U /^ U /+) are also given up to multiplications 
by phase factors, which corresponds to the action of the group U{1). The transport laws as well as 
further details of the gauge conditions are encoded in the form of the data and certain properties 
of the unknowns for the reduced equations. 

Since it turns out to be most convenient, we will carry out all our calculations on the manifold 
Ma and use for the subsets of M^ introduced above the same names as for their images under tt. 

We denote by Zui the vector field generated by Ui and the obvious action of SU{2) on M^ 
and define complex vector fields X+ = ^(^U2 + J^-ui), X^ = — (^-u^ ~ iZ^i), X — — 2iZ„3 which 
satisfy the commutation relations 

[X, X+]^2X+, [X, X-]^~2X-, [X+, X-]^-X. (2.2) 

The conformal field equations, in the form used in |^, are given in a particular (conformal, 
coordinate, and frame) gauge which is explained, together with the equations, most naturally in 
the context of normal conformal Cartan connections (cf. y]). Again, we shall not go through 
the complete argument but just describe the unknowns and equations. To obtain the equations 
on M'a, we use the fact that the solder and the connection forms on the bundle of spin frames 
induce corresponding forms (t°° , cj" f, on M^ \ /' which extend smoothly to M^. The metric 
Eab Ca'h' o""" a^^ ou M'^ is degenerate because < a"-'^ , AT > = (the angle brackets denoting the 
dual pairing), but it descends to the Lorentz metric g on 7r(M^ \I')- 

The equations are written as equations for the "vector'-valued unknown 

U — {c ab, C ab, C ab, X(ab)cd, ^abcd, fab, 0(afc)cd, ©/ah' 't'abcd), 

whose components have the following meaning. We consider the smooth vector fields 

^aa' ^ aa' ^r ' ^ aa'^p i ^ aa' ^+ i C aa' ^ — 

which satisfy < cr'^'^ , Cbb' > = Cfc " £6' ° on Af^ \ /'. All fields are written in space spinor notation 
based on the vector field V2 9r = r"" Caa' ■ Since t°''^ Caa' is invariant under the action of C/(l) it 
descends to a vector field on 7r(M^\/') which is time-like, has norm Taa' '''"■"' = 2, and is orthogonal 
to S. We have 

^aa' /= ^aa' ^r ^ a' ^ab \^'^) 

with Cab = T(a ^ Cb)b' — c° ab c'r + c^ abdp + c+ abX+ + c~ abX- . The conucction dcfiucs connection 
coefficients Tabcd — Tb"" Taa'cd — Tb"" < ujcd, Caa' > which Can be decomposed in the form 

-^ abed 1^ \sabcd \abcd} 1^ \<,abcd /(.(ab)cd} T\ ^ab Jed, 
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with fields satisfying Xabcd = Xab(cd), iabcd = i{ab)(cd)-, fab = f(ab)- The curvature is represented 
by the rescaled conformal Weyl spinor field ipabcd = (f'iabcd) and by a spinor field Qabcd — 0a6(cd) 
which is the Ricci spinor field of a certain Weyl connection for g. 

The pull back n* 8, again referred to as the conformal factor and denoted by Q, extends 
smoothly to M^ and is known in our gauge explicitly. It is given by 

and appears, together with the 1-form 

UXab-pDabU-p'^DabW 



2p- 



(U + pW)^ 



(with Xab as given in appendix |A^|) which characterizes in a certain way the difference between 
the Levi-Civita connection of g and the Weyl connection referred to above, as coefficient in the 
conformal field equations. We have set here 

" ^ (U+pW)'^' 

uj=2n{-Dab^D''''n)-^ ^ p{U + pW){U^ + 2pUx''^DabU - p^W'^UDabU (2.5) 

+2 p^ U x''^ DabW - 2 p^ D'^^U DabW - p* W'^W DabWy^ , 

where the smooth functions U = C/(p, i), W = W{p, t), which satisfy U — 1 and W = \mADM on 
/°, are given as part of the initial data on the initial hypersurface C", on which Dab is the intrinsic 
covariant derivative. Note that the fields Q., w, dab do not depend on r. The conformal factor 
satisfies the relations (cf. |Q|) 

e>o on M'a, {e = o} = y-u/'-u/'u/'+uy+, 

caa'(e)^0 and e"''e'^''''caa'(e)cfc6,(e) = on j'^. 

In the following we shall refer to the coordinates t, p, i, the frame {caa'}, and the conformal gauge 
defined by ( |2.4D as the F-gauge. 

2.1.1 The conformal evolution equations 

We recall here a few general features of the conformal field equations and refer again to [Q for 
more details. The conformal field equations imply on M^ evolution equations of the form 

{A°dr + A^dp + A+X+ + A-X^}u = Cu, (2.7) 

where A'^ , A^, A^, C denote matrix- valued functions which depend on u and the coordinates. The 
system is, for u close to the data given below and for the coordinates taking values on M^ near 
C", symmetric hyperbolic. Writing u — {v,4>) with 

V = {c ab, C ab, C ab, X(ab)cd, iabcd, fab, Q {ab)cd, ^g^ab)^ ^ = {(f'abcd), (2.8) 

the evolution equations for v are obtained, with our assumptions on the gauge, from the structural 
equations of the normal conformal Cartan connection associated with g. They read explicitly 

drc\b^-X(ab)'^c\f-fab, (2.9) 

9rC"a6 = -X(a6)'^c"e/, «=!,+,-, (2.10) 



dr^abcd = -X(ab) " £,efcd + —7= {(-ac X{bd)ef + <^bd X(ac)ef) / (2-11) 
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-\/2 X(ab)(c *" fd)e - 2 (^ac Q f'' bd + £bd 6/ ''' ac) - i Q t^^abcd, 

drfab = -X{ab)'^ fef + ^QfKb, (2.12) 

dTX{ab)cd = -X(ab) " Xefcd " 0(cd)ab + © Vabcd, (2-13) 

5r6(ab)crf = -X(cd) ''■'' 'S>{ab)ef - drO Tjabcd + l V^ (f (afJ-b)cde, (2-14) 

drOg ' ab = -X(afc) "-^ Og ^ e/ + \/2 d"-^ ry.be/ , (2.15) 

where IJabcd = \ {4>abcd + 'i'tbcd) ^^d ^afccd = -f ('/'a&cd " '/'afccd)' '^ith Ta"" U^ T^" Td"^ (jia'b'dd' 

— 0^fccd' denote the electric and the magnetic part of (jiabcd respectively. These equations are of 
the form 

drV^K{v) + Q{v,v)+L{(t)), (2.16) 

with a linear function K and a quadratic function Q of w, both with constant coefficients, and a 
linear function L oi 4> with coefficients which depend on the coordinates. We have L = on /'. 
The evolution equations for (j), derived from the Bianchi identities, are genuine partial differential 
equations. They will be considered in more detail below. 

2.1.2 The initial data 

Consequences of the finite regular initial value problem have been worked out so far for Cauchy 
data which are time symmetric and admit a smooth extension through space-like infinity. In fact, 
it has been assumed in [y, as will be done in the following, that the conformal structure is analytic 
near space-like infinity. We note that this condition is imposed only for convenience and could be 
relaxed. The free Cauchy data on S are then given by the conformal structure of a smooth metric 
h on S which is analytic in some /i— normal coordinates near i. 

We assume h to be given near z in a certain conformal gauge, the cn-gauge (cf. |Q|). This reduces 
the freedom of performing conformal rescalings h ^ 0"^ hto the choice of the 4 real parameters 9(i), 
0^a{i), the value of 6* in a neighborhood of i then being determined by the conformal gauge. We 
assume that Ba is a convex ft,— normal neighborhood of i and that p descends to a radial normal 
coordinate on Ba- 

The metric h induced by g on 5 is related to /i by a rescaling h — fi^^ h, where the conformal 
factor fl satisfies p n~ 2 ^ 1 as p ^ and the Lichnerowicz (Yamabe) equation 

(I?„I?"-ir)(r!-^) = 0. (2.17) 



Here D denotes the covariant derivative and r the Ricci scalar of h. The form ( p.SD of Q in terms 
of the functions U and M^ is a consequence of this equation and the required asymptotic behavior 
of fl, which ensures that h is asymptotically fiat. 

The initial data on C" for the conformal field equations are derived from h and fl. They are 
given by 

C°afc = 0, C^ab = pXab, C+ ab = Zab + P C^ ab, C~ ab = Vab + P C~ ab, 

X{ab)cd = 0, ^abcd^ V^pjabcd, fab = Xab, (2.18) 

&abcd ~ ^^ D(ab D cd)^ ^ vl P^ ''^ ^abcd, 4labcd — q5- {D(^ab Dcd)^ + ^ Sabcd), 
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with Xab, Vab, Zgb, and the expression habcd of the metric h in space spinor notation as given in 
appendix |A.2|, and Sabcd = S(a6cd) the trace free part of the Ricci tensor of h. 



In chapter |4J.| we shall discuss how the coefHcients c^ ab, labcd defining the frame and the 
connection coefHcients are determined on C" by the (3— dimensional) structure equations from r 
and Sabcd- The observation (cf. |Q]) that the data above extend smoothly to / ° C C" is most 
important for our construction. 

2.1.3 The transport equations on / 

At first sight it may appear that the initial data on S, thus in particular on C", should be com- 



plemented by boundary data on /' for the solutions of equations (2.7) to be uniquely determined. 
However, it turns out that for any smooth solution to the evolution equations on M'^ which coincides 
on C" with the initial data above, we have the important relation 

^^ = on /'. (2.19) 



As a consequence, equations (2^) reduce to a symmetric hyperbolic system of the form {A^ dr + 
A+ X+ + A" X- } u = C u on I' which allows us to determine the unknown u on /' uniquely in 
terms of the value of u on / ". Thus we find, as was to be expected, that any smooth solution of 



( |2.7D on M'a taking on C" our initial data is determined uniquely by its data on S. 

More generally, by applying repeatedly the derivative operator dp to the evolution equations. 



restricting to /', and observing ( 2.19 ), we obtain symmetric hyperbolic transport equations 



{A°dr+A+X++A-X^}uP = Cj,uP+gp on /', p = 0,l,2,..., (2.20) 

for the quantities u^ — {d^ u)\ji. Here the matrix- valued function Cp and the vector valued function 
gp depend on p and the quantities u°, . . . , u^~^, but the matrices A", A^ are universal in the sense 
that they depend neither on p nor on the initial data. We shall employ the notation above more 
generally, such that applying it to the fields Sabcd and r on the Cauchy hypersurface we have 
^Ibcd = idpSabcd)\j'o and rP = (9^r)|^'o, r espe ctively. 



To integrate the transport equations (2.2C) on /', we expand all fields in terms of the matrix 
elements of unitary representations of SU{2) which are given, in terms of the matrix elements 
{f^ b)a,b=o,i of the 2— dimensional standard representation of t g SU{2), by the complex-valued 
functions 

SU{2,C) 3 t ^ Tra^ ,{t) = (7)^ (7)^ t(^^ (,, . . .t^").- ,„)„ To%{t) = 1, ^^_^^^ 

J, A: = 0, . . . , m, m = 1, 2, 3, . . . 

Here, as in the following, setting a string of indices into brackets with a lower index k is meant 
to indicate that the indices are symmetrized and then k of them are set equal to 1 while the 
remaining ones are set equal to 0. The functions ^/m + 1 Tm ■' k{t) form a complete orthonormal 
set in the Hilbert space L^(/i, S'C/(2)) where /i denotes the normalized Haar-measure on SU{2). 
Under complex conjugation we have 



T™ 1 kit) = (-1)^+'- T™ ™-^" ™_fe(i), t e SU{2), 
and, for < /c, j < to, m = 0, 1, 2, . . . , we have with pm.j = {j (to, — j + 1)} 2 

A function / satisfying a relation Xf — 2sf with an integer or half integer number s, is said to 
have spin weight s. We note the spin raising (lowering) property of the action of X± on such 



functions implied by (2.2), i.e. X X± / = 2 (s ± 1) X± f. By construction of the manifold M„ any 
function occuring in our formalism has a well defined spin weight. This leads to a simplification 
of the expansion in terms of the functions T^^j . The general form of these expansions has been 
discussed in detail in 0] and will be assumed here without further explanation. 
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The quantities m°, u^, u^ have been determined in [Q. They are given here (with a correction 



and a useful change of notation) at the beginning of chapter |0|. The functions u^ will be 



calculated in chapter |4J.|. The quantities u^, p = 2, 3, . . . have been shown (cf. ^) to develop a 
certain type of logarithmic singularity on the sets / * unless the free datum h on S satisfies the 
asymptotic regularity condition 

Dia,b,---Da,Hbabcd){i)^0, (2.23) 

for g = 0, 1, 2, ... , where the spinor field babcd = b(^abcd} represents the Cotton tensor of h. The 
values of the func tions u^, p < 3, which will be given below, have been calculated on /' under the 
assumption that ( 2.23| ) is satisfied for q < 1. The analysis of the quantities u^, to the extent to 



which it has been carried out in y|, indicates another potential sour ce fo r a singular behavior of 
the fields u^, p > 3, at / *. This will be discussed further in chapter |4.1|. 



2.2 The NP-gauge 

For simplicity we restrict our discussions now to the future of S* in M, we refer to future null infinity 
simply as to null infinity and we denote it by J. In the following we shall describe a certain class of 
gauge conditions on (Af , g) near null infinity, referred to as the NP-gauge, which comprise certain 
requirements on the conformal gauge, certain coordinates, and a certain orthonormal frame field. 
Though this gauge is known, our description will be quite detailed, because we will have to refer 
to it later. The Levi-Civita connection induced by the conformal metric g will be denoted by V. 

Suppose {E°^,} is a smooth frame field, satisfying g(E'^a' ^ -^w) ~ ^ab^a'b', which is defined 
in a neighborhood of null infinity. We call it an "adapted frame", if it satisfies the following 
conditions. The vector field E^^, is tangent to and parallel propagated along null infinity. On the 
neighborhood on which the frame is given there is exists a smooth function u° which induces an 
affine parameter on the null generators of J such that E°j^, {u°) = 1, which is constant on null 
hypersurfaces transverse to J7, and which satisfies Eq^ — g^^S/ pu° . Thus Eqq, is tangent to the 
hypersurfaces {u° = const.} and geodesic. The fields E°j^,, Eqq, as well as the fields i?oi'' -^lo'' 
which are necessarily tangent to the slices {u° — const.} n J, are parallelly propagated in the 
direction of E^^, . 

In terms of its NP-spin-coefHcients (note the slight difference of our notation with that of ||ll|| ) 

^la'bc = 2 \Ka'Kl'^aE°Q,f, + E°^,E°yVaE^Q,i^j , (2.24) 

an adapted frame is characterized by the properties 

rVii = o> rji,ii = o on J, 

^ lo'oo ^ ^ oi'o'o'' ^ ii'oo ^ ^ oi'o'i' + ^ oi'oi' ^oo'a6~"' a,,0 — v,l, near j. 



(2.25) 



The first of these conditions tells us that J^ is shear free. This well known fact follows from the 
equation for the trace free part Sa(3 of the Ricci tensor of the conformal vacuum metric g, 



e s^p = :r5a/3V^vTe - 2v„v^e. (2.26) 



Transvection with E^f^E^^^ and restriction to J gives Tl^.n E^g,{e) = 0, while -Bgo'C©) ^OonJ. 
We shall combine now the construction of an adapted frame with the freedom to perform rescalings 



q2 



g, e~>e* = 6ie (2.27) 

with some positive func tion 9, to obtain another adapted frame {E'^,} for which we get further 
simplifications besides ( 2.25| ). We start with an adapted frame {E°^,} as described above. For 



arbitrary 9 > and for arbitrary function p > which is constant on the generators of J we set 
Ely ^ 9-"^ p El^, and u'{u°)^ 9^{u')p-^{u') du' + ul on J, (2.28) 
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where the integration is performed along the generators of J'. Then E'y will be parallelly propa- 
gated and E'y{u') = 1 will hold. We assume that u° = u° and u' — ul on C and set 

E'o„,^p-^E°„,, El,,^9-^pE°,„ E'o,,^0-^E°,, on C. (2.29) 

Since C is difleomorphic to S'^ and thus carries (up to diffeomorphisms) precisely one Rie- 
mannian conformal structure, we can fix coordinates x^ = "d, x'^ — (p as well as the func- 
tion 6* on C such that the metric h* induced by g* on C is given by the standard 5^— metric 
h* = dd'^ + siia^ d d ip'^ . Using the transformation laws r*Q,QQ — p~^ [rio'oo ~ -^oo'(l°S^)] ^^^ 
-'^oi'ii ~ P&^^ [rQi'ii + -^11' (log ^)] oil ^1 we can achieve, by suitable choice of d^ and p, 

rVoo-0, r*i,ii = 0, £;oV(0*) = const. ^0 on C. (2.30) 

The transformation s*^ = -| {(VaV/jS - ^Vad^pd) - hapj^-y^y - lV.^9V^e)] + s^p of the 



trace free part Sap of the Ricci tensor under the rescaling ( |2.27D implies a transformation of 



<&22 = \sapEly Ely into $22 = 5*Q/3-^ri"^ii^i which yields, with the assumption that $22 = 
on J, on the generators of J the ODE 

Ely [Ely {9)) - ^(i?°i,(0))'-e$22 -0. (2.31) 

This equation can be rewritten as a linear ODE for 9~^ which can be solved on the generators of 
J with 9 > Q. Using the initial data 9, Ely {9) on C determined above, we solve for 6 to obtain 

*22 = 0, r*i,ii = on J. (2.32) 

Here the second equa tion is a consequence of the first, the field equations, and ( ^.30| ). We assume 



in the following ( 2.28 ). We observe that the induced metric on the sections {u* = const.} is given 
as a consequence everywhere on J^ by the S'^-standard metric. 

Once 9 and E^y have been fixed on J, the vector field E'y (whence E^q,) tangent to {u* = 
const.} is determined up to rotations. We choose some smooth field E'y on J, solve the equation 

Ely (c) = -t El,^ Elf, V^ E'oy^ (2.33) 

for the function c with initial value c = on C, and replace E'y by e^'^ E'y to achieve 

rti,oi = on J. (2.34) 



Observing the simplifications above, we contract the analogue of ( 2.26 ) for g* with E'yEl^, to 



conclude that V^V^^O"^ = on J^. A further contraction with E'ffEly gives 

Elr{E^o'i^*)) "^^^ i-e- E^^'i^*) ^ const, on J, (2.35) 

while a contraction with E'^E'y yields now E'y (^E'q, (Q*)^ ~ TlyQQE'Q,{Q*), which implies 

r*i,oo = on J. (2.36) 

To fix also d9 on J, we use the conformal transformation law for the Ricci scalar, i.e. 

i?[5*] = li?[.g] + ^V'Lf?V*"0-^V*„V'^"0. (2.37) 

If we require that R[g*] — along J, this equation takes on the generators of the null hypersurface 
J the form 

Elv {Ko' (0)) - hly i9)E',o, (9) = F\ (2.38) 
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of a linear ODE for the unknown E'q,{9), where the right hand side 

F* = 9\c i^E'oy{E',o,{e)) ~2T',,,,,El,,{e) - ^E'oy{e)E',o,{e) + ^R[g] 
is given in terms of quantities which have been determined already on J'. Using the initial value 



-^00' (^) — P ^r°o'oolc, fixed on C by ( |2.3(]| ), we can integrate the equation to achieve 



R[g*]^0., r'o-oo-O on J, (2.39) 

where the second equation follows again from our previous results and the field equations. 

We do not require conditions of higher order on the conformal gauge. Assuming a conformal 
gauge as described here, we shall refer to an adapted frame {E'^,} satisfying the conditions above 
as to an NP-frame, and to a normalized spin frame e*'^ = {o'^,l'^} which implies a NP-frame 
as to a NP- spin- frame. 

We extend the coordinates x^, x^ to J such that they are constant on the null generators 
of J . As described above, we define null hypersurfaces {u* = const. ^ transverse to J and we 
denote by r' the affine parameter on the null generators of these hypersurfaces which satisfies 
E^Q,(r*') — 1 and, on J , r* = 0. The coordinates x^ , x'* are extended such that they are constant 
on the null generators of {u* — const.}. Thus we get a Bondi-type system (w*, r*,x^, a;'') in some 
neighborhood of null infinity. Occasionally we shall change from the coordinates ■(?, Lp, to a complex 
stereo-graphical coordinate given by C = e'^^'^ctg^. We write the volume element and the volume 
form alternatively 

ds^ = -(rfi?2 + sin^^dip^) = -P{Cy^dCdC, e = sin^d^ A d(p = [2F(C)]"^<iC A dC, 

where we set P{() = ^(1 + CO- We shall refer to the conditions on the conformal scaling, the 
frame field, and the coordinates as to the NP-gauge. 

2.3 Relating the NP-gauge to the F-gauge 

While the NP-gauge is hinged on null infinity, the F-gauge is based on a Cauchy hypersurface 
and these gauge conditions are in general completely different. In the following we will study the 
transformation which relates one to the other. It is important for this that the conformal factor 
6, whence J, is known explicitly in the F-gauge. 

The vector fields {caa'} tangent to the 5-dimensional bundle space M'^ are not directly related 
to the NP-gauge on the subset Mo \ / of M . Let 5^ D C/ 9 p ^ s{p) e SU{2) be a smooth local 
section, defined on some open subset U of S*^, of the Hopf fibration SU{2) -^ SU(2)/U{1) ~ 5^. 

It induces a smooth section C/ x M x M 9 {p,T,p) -^ {s{p),T,p) £ M'^. We denote the image of 
S by M*. The vector fields tangent to s{U) which have projection identical to that of X± are of 
the form X± -\- a±X with some smooth functions a± on s(f/), satisfying a_ — —a+. Because of 
( |2.2D a± cannot vanish on open subsets of s{U). Consequently, the tangent vector fields c*^/ of 
M* satisfying tt* (c*„, ) — tt* [caa' ) are given on M* by 

'^*aa' = Caa' + («+ C+ aa' + «- C" aa') X , 

with functions a± which are independent of r and p. The connection coefficients defined on M* 
by the connection form oJ' c and the vector fields c*^/ are given by 

r:,, ' c = r,,, \+{a+ C+ aa- + a- C- aa') (eq ' Ec " " e/ Ec ')• 

In the remaining part of this section we shall work on n{M'a) and denote the projection of the 
vector fields c*^,,, which define a smooth orthonormal frame field on 7r(M* \ /'), and the pull-back 
of r*|j, ^ c by S" again by c*^, and r*^^, ^ c- Similarly, the projections of J^' n M* and / + n AI* will 
be denoted by J^ and /+. 
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The frame field {c*^/}, which is in general not adapted to null infinity, will now be related close 
to /■'" to an adapted frame {E^^'}- On J' the vector field E°y must be of the form 



E°{^ = /v"e, 



(2.40) 



where V and 8 denote the Levi-Civita connection and the conformal factor associated with the 
F-gauge. The requirement = El^,V pE^p = fV^QW^f V^G + pV p {\ V„e V"e) that E^p 
be parallelly propagated, gives after contraction with a vector field Z transverse to J the ODE 



v"eVa(iog/) = - 



zW) 



(2.41) 



for / on the generators of J. To fix /, we set f = fo = const. > on some section C of J'. The 
function u° satisfying E^{f{u°) = 1 on i/ and u° — u° on C can be now be determined. 
Let A'^f, e SL{2, C) satisfy 



rpo \b \b' * 

^aa' — ^ a^ a' ^bb' 



(2.42) 



Rewriting ( |2.4(]| ) in the form E°y = /c^j,(9)e°''e'''^'c*a,, we find the relations 

x\x°\, = /c*i,(e), x\x^\, = -/c*o,(e), x\x^\, = fc*,,,{e). (2.43) 

From ( |2.42| ) we obtain X°iE°j^, = A'^g-^'ii' ~ ^"'I'^io' ~ A^'-^,c*;^,. Applying this to the function u° , 



we get 



x'o^~^\'Cio,{un + x\,ci,,{un. 



(2.44) 



Together with the condition det(A"j,) = 1 the relations (2.43), (2.44) allow us to determine the 
matrix elements X"^^ on J up to replacements A°{, —t A"j r]"^ with [rf b) = diag{e^",e~^") G U{1). 
After making here an arbitrary choice, the adapted frame {E°^,} is determined uniquel y nea r J. 
To determine an NP-frame {E*^,} near J, we need to find an appropriate rescaling ( 2.27 ) and 
a scaling factor p. We set 






E', 



A'' A*-' 



a' ^bb' 



with A"t e SL{2,C). Assuming (2^, we have E'^ = /* V*"e* with 



r = ^ and E',,,ie*) = j^ on J. 



(2.45) 



(2.46) 



We choose now 9, dO, and coordinates x^, x^ such that the induced met ric on C is given by the 
S'^-standard metric and, with -p chosen such that p = on C, conditions ( 2.30| ) are satisfied with 

Following the p rocedure of the previous section, we can determine the conformal factor Q on 
J such that ( 2.32| ) is satisfied. The transformation A^^ can be determined in the same way as 
A"j. Imposing condition (2.34), we determine A'^j, up to C/(l)— transformations on C. Conditions 
(2.35), (2.36) will now be satisfied as well and we can determine dO on J such that (2.39) holds. 
Extending the tetrad to a neighborhood of J such that it is parallelly propagated in the direction 
of i?oo'' ^^ S^t the desired NP-frame. 

In our later calculations we will need the quantities i?QQ, (A°^). Using our gauge condition 
^oo'afc ~ and the transformation laws for the connection coefficients. 



F* 






ea(6C, 



)a'(l0ge)}, 



El^. (A" ,) = -A/ a A^' a' A" , F* , \ + A^ F^, 
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where r*^,^,^ denotes the connection coefficients with respect to V* and {c*^/}, we find 

E'o„, (A" ,) = -A^ A^' 0' A'' , r}^, " ,. (2.47) 

In the considerations above we had to fix various quantities by prescribing data on the section 
C. When we shall determine later the expansion of a NP-frame near /+, it will be natural to try 
pushing C to /■'". A priori it is not clear, however, whether this can be done in a continuous way. 
We shall see, that for certain quantities the limits to /"'' do exist, while others quantities can only 
described in terms of their growth behavior near /+ . 

3 The NP-constants 

In 1965 Newman and Penrose discovered certain non-trivial quantities, defined by certain integrals 
over a 2-dimensional cross-section of J'^ , which are absolu tely conserved in the sense that their 
values do not depend on the choice of the section (cf. ||l2|, |l4l). The interpretation of these ten 
real NP-constants is still open. In the case where the space-time admits a smooth conformal 
extension containing a point z+ ("future time-like infinity") whose past light cone represents J^ , 
these constants are essentially given by the five complex components of the rescaled conformal 
Weyl spinor (cf. ||l4| , 0|). However, these quantities do not allow us a simple interpretation either. 
More interesting is the case of stationary vacuum space-times. In this case the constants have been 
calculated and have been given in the form {mass) x {quadrupole moment) — {dipole m,oment)^ 
(cf. 0, @). 

If the evolution of the field in time is not given explicitly as in the presence of a time-like 
KiUing vector field, there appears to be no obvious way to calculate the NP-constants. It turns 
out, however, that under suitable assumptions on the asymptotic behavior of the field near space- 
like infinity the constants can be calculated by integrating the transport equations on /' to a 
sufficiently high order. In the following we shall derive a formula for the constants in terms of 
quantities which can be determined by solving the transport equations. 

To explain the original formula (cf. IQ), which is given in the Bondi-Sachs-Newman-Penrose 
framework, let (u, r, ^, ip) denote Bondi-coordinates on the physical space-time, where r denotes an 
affine parameter along the generators of the null hypersurfaces {u = const.} and the generators are 
labeled by the standard coordinates (^, v?) on the two-sphere. The null frame {Eaa'} as well as a 
corresponding spinor dyad {5"^, l^}, both defined on the physical space-time, are normalized with 
respect to the physical metric g. They are adapted to the Bondi-coordinates such that Eqo' = 9r- 

We assume that the conformal space-time with metric g* :— r~^g admits a smooth extension as 
r ^ oo to a smooth Lorentz space with boundary J'^ = {r* = 0} and that the functions u' :— u, 
r' := r^^, ■&, and (f extend such as to define a smooth system of Bondi-type coordinates near »/+. 
Furthermore, we assume that the frame {E'^^,} and the spinor dyad {o'^, t*^}, defined by 



E'a' 



^2 — a— a' 

3.1) 



such that they are normaHzed with respect to g*, extend to smooth frame resp. dyad near »/+. 
The results of Newman and Unti (cf. ||l5|) then imply that {E'^,} defines in fact a NP-frame. 

Under our assumptions the component -00 = iPabcdo'^o^ o'^' o^ of the conformal Weyl spinor 
has an expansion V'o = V'o''"^ + V'o''"^ + 0{r~'^) with coefficients V'o which are independent of r. 
In terms of the physical space-time the NP-constants are given with this notation by the integrals 

Gm = f 2Y2.,m'lpl sIm? M dip, (3.2) 



which are calculated for fixed value of u. The functions 2^2, m, Tn — —2,-1,0, 1,2, denote spin-2 
spherical harmonics (cf. |^) which are obtained from the standard spherical harmonics by 

2>^2,m = — --=£'oi"-Boi''5a'^/3^2,m = yif^ ^2,m- (3.3) 
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Here S and 5 denote the standard covariant differential operator on the unit 2-sphere and the 



"edth"-operator, respectively. In evaluating (|3.2D , it will be important that the operator 8 is 



defined with respect to the complex null vector field E'y (cf. |10|| ). 

We reex press the constants in terms of the fields g* , E'^^,, o , l'"^ satisfying the NP-gauge, in 
particular ( |2.39| ). Using the component 0o ~ r iI^abcdo' ^o" ^ o" ^ o" ^ of the rescaled conformal 
Weyl spinor, and performing the obvious lift to M\ we obtain for the NP-constants the formula 

G„ = -^ j>2Y2..mE'Q^,{(t)o)dSda. (3.4) 

Here dS ~ smd d-d dip denotes the surface element on the cross-section {r*, u' = const.} C J^ and 
a denotes a parameter on the fibers of the principal fiber bundle M' -^ M . The second integration 
can be performed without changing the result because the integrand is independent of the variable 
a. 

The values of these integrals are independent of the value of the constant defining the cross- 
section as well as of the choice of the Bondi-coordinate u' itself. Thus they are invariant under 
supertranslations (cf. |Q). 

We shall determine the NP-constants by integrating the transport equations on /'. Since these 
equations and their unknowns are given in the F-gauge, we express (6^4) in this gauge. Using 



( 2.45 ), we obtain in the notation of the previous chapter 

G™ = -^^2i'2,™^{A%A%A%A%[A%AVcL,(^,,rf,)-3 0,,,,i^oVW] 

(3.5) 

+ AeK\K\K\E'^^,{K\)4>b,de}dSda. 
This is the expression for the NP-constants which will be used in the calculations of section |4 



4 Time symmetric space-times 

In this section we will use the assumptions of the regular finite initial value problem near space-like 
infinity and thus restrict our considerations to time symmetric space-times. We begin by solving 
the third order transport equations on /'. This calculation is of interest for two quite different 
reasons. First of all, it will give us a first insight into the potential source of singular behavior of 



the quantities u^ pointed out in section |2.1.3|1 . Further, besides giving information on this question 
of principle, the calculation will allow us to analyze the relation between the NP-constants and the 
initial data for asymptotically fiat solutions. Under our assumptions, we will be able to evaluate 
the integral (^^) in terms of quantities derived from the initial data. 

4.1 Solving the third-order transport equation 



The solutions u^ of equations ( p. 201 ) have been given in ||4[] for p < 2. Since they will be used in the 
following calculations we reproduce them here, in a notation, though, which is more convenient for 
a systematic discussion of the higher order expansion coefficients. We also take the opportunity to 
correct a misprint in 



The solution vP of the transport equations (2.20) has the form 



where m — niADM denotes the ADM-mass of the initial data set. The spinors appearing on the 



right hand side of these and the following formulae are listed in (A.IC) of appendix |A.3|. The 
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solution u is given by 

(C^b)^ = C±l(r)yafc, $1^^^ = 5'HT)(eaca;fcd + eMXac), Xlab)cd = ^HT)£abcd' 



fL^F'ir)xa,, iQ',)'^tHr)xa,, 



e 



{ab)cd 



THr)el 



(4.2) 



bed' 



while u^ takes the form 

ic"abf = [cf (r) + cf(T)Wi]Xab + cf {T)[X^Wiyab + X+WiZab], 

iC'\bf =C^'^{T)Xab, 

i<bf = [cf W + 4^{r)W^]Zab + 4^{T)X_WiXab, 

Kb)^ = [4'ir) + 4Hr)W,]yab + 4HT)X+WiXab, 

Cbcd = [5'?(t) + S2{T)Wi]{eacXbd + e.bdXac) + Sl{T){eacybd + £bdyac)X^Wi 

+ Si{T){eacZbd + ebdZac)X+Wi + SliT){el,,aX+Wi + sl^.^X-Wi), 

Xtab)cd - [Ki{t) + K^iT)Wi]sl,,a + Kl{T)habcd + Kl{T){eacybd + ebdyac)X^Wi 

- Kl{T){eacZbd + ebdZac)X+Wi + Kl{T){el^^^X+Wi - e^.^^X-VFi), 
/ab = [Fi{r) + F|(T)W^i]a;,b + F^{T){X^Wiyab + X+T^i^ab), 

(0/ab)' = [<? W + ^2 WM^ll^ab + t2(r)(X_W^i2;,fc + X+W,Zab), 

[Ti{r) + TiiT)Wi]sl,,a + T^{r)habcd + Tl{T){eacybd + ebdyac)X.Wi 
" Tl{T){eacZbd + ebd;^ac)X+W^i + Ti{T){e\,,^X+Wi ~ elf,,^X^W^), 



(4.3) 



«2 



-/-abed = </'?(T)X+X+W^2eLd + [-^2(^)^+1^1 + 4>1{T)X+W2]e\,,a 

+ [cfiir) + ct>l{T)W^ + 0i(r)M/2]eLd " [4>l{-r)X^Wr + 02(_r)X_Ty2]eLd 
+ </.2(-T)X_X_T4^2eLd- 



The r-dependent functions in these expressions are polynomials which are given in appendix! A. 3|. 
The calculation of u^ is facilitated by the following properties of the transport equations ( 2.201 ) . 
For p>l they are of the form 



drV^ ^LpVP+ Ip, B'^da(l)P = Mp(j)P, 



(4.4) 



where, using the notation (2^), we set v^ = (S^w)!/', 0^ ~ {dP(l>)\i> and denote by Lp and Ip a 
matrix- resp. vector-valued function of the quantities u°, . . . ,uP~^, while Mp denotes a matrix- 
valued function which depends on the variables u", . . . , u^^^, v^. The matrices B°' neither depend 
on p nor on the initial data. Thus, given the quantities u', g < p — 1, we can integrate the first of 
equations (4^), which is an ODE. To i ntegr ate the sec ond equation, we expand the quantities u^ 



in terms of the functions T„ ^ given in ( 2.21 ) and use ( 2.22 ) to reduce the integration to that of a 
system of ODE's. 



To determine the initial data for u^ on I ^, we have to expand the unknowns ( 2.18 ) in terms 
of p. Instead of prescribing the conformal metric h on the initial sHce, which represents the free 
datum, we shall prescribe, in a fashion consistent with the 3-dimensional Bianchi identities, certain 
curvature quantities and use the 3-dimensional structure equations and the Yamabe equation to 
determine the remaining quantities. 



The conformal factor, which appears in the expressions ( p.lq ), is given in (2.5) in terms of the 
functions U and W. The function U, which is determined locally by h near space-Hke infinity, is 
given, by a procedure explained in f^, in the form 



oo 



p=0 



(4.5) 
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with p-dependent coefficients Up. As shown in ^, the Taylor expansion of U in terms of p has in 
our gauge the form 

For our calculations we shall need the coefficient U4, which will be determined later in this chapter. 
The function W, which contains global information on the free initial data, is determined by 
solving the Yamabe equation on the initial hypersurface. We shall consider here a larger class of 
functions which are subject to the Yamabe equation only in a small neighborhood of space-like 
infinity. The coefficients in the Taylor expansion W = Wq + Wi p+^W2 p^ + ^Ws p^ + 0{p^) have 
expansion (cf. |Q|) 

2i rn 



m=0 fc=0 



W^^}^}_^ W^,^.k T, 



They are restricted by the requirement that the Yamabe equation {h^^DaDfi — ^rii)[W] — holds 
near {p = 0} , which implies the simplification 

2 
W,=Y,W^■,2^,kT2^, z<3. (4.7) 

fc=0 

We get for the conformal factor and the trace-free part of its second covariant derivative 

n^p^ ^mp^+ [| m2 - 2 Vl^i] p4 + [-^m^ + 3mWi - W2] p^ 

D^abD^d)^ = [~6mel,Jp+ [(12 m^ - 36W,)sl,^, 12 iel,^,X+ ~ sl,^,X^)W,]p^ 
+ [(-15m3 + 96mW^- 36W2)el,^^i + {ei,,aX+ - e3,,,X_)(24m W^i - 8W2) 
-U^lbcd^+^+ + ^lbcd^-X-W2]p^ + [{156 Wi - 150 m^Wi + 15 m^ + 81 mW2 (4-8) 
-20 W3 -AU4 + ^X+X^Ui -6X+WiX^Wi)el,^a 
+ ieLcdX+ - eLd^-)(30M/2 - SOm^l^i + 15mW2- fW^ - lU^) 
+ WlbcdX+X++st,^,X^X.){SW? + lmW2- iW^3 - i^f/4) - h,ial\,^,)]p' 
+0{p'). 

From this we obtain as initial data for m'^ on / '^ 

(O'=o, {c\,r^o, ic\,f^o, ic-,r=o, 

^Ucd= 3^+^+W^2eLd + (-72 mX+W^ + i8X+W2)sl,^, 
+ (27 m3 - 288 mWi + 216 T^2)eafccd 

+ i72mX^W4 - i8X^W2)el,^^ + 3X^X^W2ei,^,, (4.9) 

^Ibcd = i^lbcdX+X+ + et,^,X^X^){9Wl lmW2 - W^ - ^,U,) 
+4{el,^X+ - el,,,X^)i9W^ - f mM^2 - 5W^3 - ^Ud 
+6 e^bc<i(12 Wl -3mW2- 2OW3 - 4C/4 + y^X+X_C/4 - 6X+WiX^Wi) 
-4a:e(a7L'd) + 3s^bcd> 
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where 'jabcd — {"2^ p)^^ {^acXbd + ^bd^ac) + labcd denote the connection coefficients on C . 

We determine now how the functions t/4, labcd^ ^^^ ^abcd ^^^ related to the free data on the 
initial hypersurface C". As shown in [Q, the structure equations on C", which relate the connection 
coefficients to the curvature, read 

-^{dplOQab + — [70000^ab — JOOllVab + yjTOOab] } = 700007llab — 7001l700ab — 2^ab00 — ^y| ^ yab, 
-^{dpjllab + ^[jlWOyab - lllUHab + ^7lla6] } = 7ll007llafc ^ 7llll700afc + ^Sabll - -^ r Zab, 

and the components of jabcd have Taylor expansions 

lOlab = 0, 700ah = ^700a6 P^ + 0(p'^), 7llafc = -^iflab P^ + 0{p^). 

From this we get 

-3 32 -3 32 -3 32 

70001 — ~5y2*"0"l' 7ll01 ^ :iy2*0111' 7oOOO = ^^^*0000' 



'1100 - 5V2*0011 10V2 ' '0011 - 5V2*0011 + i0\/2 ' '"" ~ .;,/^*llll' 



^,3 _ _3_ „2 
5\/2'" 



and obtain thus for the quantity Fabcd — —4 x^^alXfd) + ^ *abcd the concise expressions 



Fo^l-sl Fi=3s?, F2 = f,s2-^r2, F^^Ssl F4 = f s^, (4.10) 



where we set Fi — -F'(afccrf)i) ^^ — S(af,cd) ) using the notation introduced in ( 2.21 ). 

In the cn-gauge the curvature vanishes at zeroth and first order at space-like infinity. At second 
order this is in general not true and the prescription of the free data on S in terms of curvature 
quantities has to be consistent with the cn-gauge, the Bianchi identity, and the regularity condition 



|2.23) for (J = 1. The content of the cn-gauge is expressed in second order in the curvature by the 



conditions 



DabD''''r = 0, DabD'^^'Scdef^-^DcdDefr, D^abDcdSefgh) =0 at I. 

It follows that the spinor 

^abcd efgh — ^ab ^cd ^efgh TT llabcd ^h ^^efghj 

where A^ denotes the Laplacian corresponding to the metric h, is symmetric in the first and the 
last four indices separately. Using the Bianchi identity 

D'^^'Sabcd^lDcdT, 

6 



we thus get 



- Dab Dcd r — - A/i Sabcd — f^ ab cdef = ta*^ b cdef — Da'^ Db'' Scdef + 7 A/i Sabcd, 

6 3 6 



whence 



19 

Da" Db^ Scdef = -^ Dab Dcd r. 

No further conditions are implied at i on the Ricci scalar r at this order. Finally, we get from 



(2.23) for g= 1 



F>^ (a F)bc Sdef)h = at i. 

The relations above imply that the expansion of tabcd efgh in terms of symmetric spinors and eafc's 
can be expressed completely in terms of symmetrized twofold contractions of this spinor, which in 
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turn can all be expressed in terms of the symmetric spinor Dab Dcd r. Working out this expansion 
we get 



jjab jjcd g 



cfgh 



h'^^\,fD-^\t,)r-^h'^^''''DefDgnr at z, 



15 



(4.11) 



in our gauge. Going through the procedure described in section (3.5) of ^ we get S(^abcd 
g2 „2 _L n( „'i\ o^^ ^ — ^2 ^2 j^ n^o3^ 



s] p2 + 0(p3) and r = H p^ + 0{p-^) with 

4 



„2 _ 3 



2-j 

12 



E^Ur'^4' -'-TaE^^^' 



'4 2j 



(4.12) 



fc=0 



fe=0 



1 /4\ 2 



where we set Rl — ^ (^) ^ D(^abDcd)^ r* , with the star indicating that the quantities are given in 
our gauge at i. The 5 real numbers R^ contain precisely the information on the metric h which 
can at this order be freely specified in the cn-gauge. 
We note that the Cotton spinor is then given at i by 



Dabbcdef = -g {^aibDcdDef) T + eb(aDcdDef) r} 



and the deviation of h from conformal fiatness at i is encoded at this order in the symmetric spinor 
Dab Dcd r {i). 

From ( |4.10| ), ( t4.12| ) we obtain 



P — 27 \ ^ r>* rp k p _ 3\ ^ TD* rpk jp _ 3 \ ^ p* rp k 

-^0 - wZ^^k ^4 0' ^1 - sZ^^k ^4 1: ^2 - ^^ /^^k ^4 2: 



fc=0 
4 



fc=0 

4 



fc=0 



TTi _ 3 \ ^ p* rp k p _ 27 \ ^ T>* rp k 

-f^3 - 8Z^-"'fe^4 3: ^4 — 2oZ^^fc-'4 4- 



fe=0 



fc=0 



Finally, we will calculate the coefficient C/4 in the Taylor series (4.6). Only the coefficients Uq, 
Ui and U2 of the expansion (4.5) contribute to t/4. These functions have the following expansions 
(cf. for the defining integrals). 



Uo ~ exp 
where we used the expansion 



{i i' (^^" + ^) 7-} = 1 + i! [^^""J P' + ^(^')' (^-13) 



Ap^ = -6+^7?ioo/ + 0(A 



Further we have, with L denoting the Yamabe operator. 



2p Jq Uo 2 



I^^3 „1,2 

36 ^""° 48 



p2 + 0(p3). 



Finally, observing ( tl.l2[ ), we obtain 

C/2 = 



C/o /•" L[C/i] p' ^ , 



2p^Jo Uo 



dp' = 0{p). 



Collecting results, we arrive at the expansion 



^ = ^+4! 



_!^^3 „1,2 

o 7iioo ^ ' 



p4 + 0(p^) = l + i 



10\/6 



E^^^^ 



*:=0 



(4.14) 



p' + Oip'). (4.15) 
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Since the initial datum for the conformal Weyl spinor is a non-Hnear function of the basic 
quantities and the transport equations are quadratic in the unknowns, we have to make use of the 
Clebsch-Gordan expansions of products Hke Tr^^T2n- These are readily calculated by using the 



definition (|2.21| ). For the quantities relevant in our calculation we thus obtain 

4 4 

X- WiX+Wi = -Y^ akTi\ + 2 6, W^ = Y. ''^^i^ + ^' 

fe=0 /c=0 

4 4 

W^X^W^ = -^Y. °fe^4''3> W^X+W^ = ^Y. «fe^A, (4.16) 



fc=0 k=0 

4 4 

{X-Wif = V^ akT^\, [X+Wif = VqY. «'«^A, 

fc=0 k=0 

with coefficients 

03= 75W^l;2,2W^l;2.1, ^4 = 7fl^l';2,2, ^ = - 1 (W^l;2,0 W^l;2,2 - ^ W^l%,l). 

It was shown in M that the quantity 0f has an expansion of the form 



Using the results above in the last equation of (O), this expansion reduces to 



(4.17) 



E ECa-^m^-2+,- (4.18) 

m=\i-2i\ k=0 



4>%m,k = 0, for i = {0, . . . , 4} and m > 8, 

<^3;6,fc =-10^/3Ty3;6,fc, 0^ g ^ = -2 ^30 W^3;6.fc , 

'/'o;4.fc = 18V6 afe - 376 m P^2:4,fc + |i?fc, 0?.4,fc = 9^/6 flfc - | V6 m W^2;4,fe + |i?fe, (4.19) 

</'2;4,fe = ISflfc - 3mW^2;4,fc + ^i?fe, 0^;4,fc = ^^^k ~ |\/6 TO VF2;4,fc + |i?fc, 

•^Ufe = 18\/6afc - 3V6TOVF2;4,fc + fi?!-, 
Given these data on / °, we are in the position to solve the transport equations on /'. The first of 
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the systems (|4.4D can be integrated step by step with the result 

{cl,f = [cf (r) + cf{T)W, + cf{T)W2]Xab + [cf{T)X+Wl + cf {T)X+W2]Zab 

+ [cfiT)X^Wi + cf{T)X^W2]yah, 
ic\bf - [clHr) + cl''{T)Wi]Xab + cmT)[X+WiZab + X^WiVab], 
(C+J' = [c±3(r)X_VKi + 4''iT)X^W2]Xab + [cfHr) + 4''{t)Wi + cf"" {T)W2]Zab 

+ 4''{T)X^X^W2yab, 

Kb)' = [cf{r)X+W^ + cfiT)X+W2]xab + [cfir) + cf{T)W, + cf (r)T4^2]j/,6 

+ cf{T)X+X+W2Zab, 

ilbcd = Sl{T)X+X+W2el,,a + [SI{t)X+W^ + Sl{T)X+W2\el,,a 

+ [SI{t)X.W^ + Sl{T)X^W2]el,,a ~ ^f (r)X_X_T^2eLrf 

+ [SI{t) + SI{t)W^ + Sl{T)W2]{eacXbd + ^MXac) 
+ [S^j{t)X+Wi + Sl(T)X+W2]{eacZbd + eMZac) 
+ [S^{t)X_W^ + Sl{T)X_W2]{eacybd + eMVac), 

xU)cd = Kl{T)X+X+W2el,,d + [KI{t)X+W^ + i^|(r)X+M^2]£Ld 

+ [Klir) + KI{t)W, + Kl{T)W2]el,,,i - [KI{t)X_W, + Kl{T)X_W2]el,,a (4.20) 
+ Kl{T)X-X-W2el,,a + [K^7{r) + Kl{T)Wi]habcd 
+ [KI{t)X^W^ + Kl^[T)X^W2]{eacybd + eM2/ac) 

- [i!:|(T)X+W^i + KUT)X+W2]{eacZbd + ebdZac), 

fib - [Ff{T) + i^lMW^i + Fl{T)W2\xab + [FI{t)X_W^ + Fl{T)X_W2\yab 

+ [FI{t)X+Wi + F|(T)X+iy2]Za6, 
(0/aJ' = [tli^) + tli^Wl + tl{r)W2]Xab + [tl{T)X_W, + tl{T) X _W2]y ab 
+ [tl{T)X+Wi + tl{T)X+W2\Zab, 

Q\ab)cd = T^{r)X+X+W2e'i,,^ + [Ti{T)X+Wi + Ti{r)X+W2]el,,, 

+ [Tiir) + Ti{T)Wi + T|(r)M/2]eLd " [Ti{T)X^Wi + Tl^X^Ty^leLd 
+ T^iT)X^X^W2el,, + [T^{t) + Tiir)Wi]habcd 
+ [Ti{T)X^Wi + T^o{r)X^W2]{eacybd + ebdyac) 

- [Ti{T)X+Wl + T^o{T)X+W2]{eacZbd + tbdZac). 



The r-dependent functions in these expressions are given in appendix ||A.3|| . 

We now turn to the second of the transport equations (14), which is a partial differential 
equation. The system for the expansion coefficients (jy^ of the rescaled conformal Weyl spinor on 
/' has the form 



(1 



T)drCl>i+X, ,^ 
drCJ^l + \X_4,l 
dr(t>l + hX^(t)l 



-(1)1= i?, 



2^4 



-2X+<i>l 



3 -^ 2^- '^2 

3 I V ^3 I j,3 



\x+4>l 



^1 = 
R21 



i?i, 



^3: 



(4.21) 



(l-T)9,0t + X_0i + 03= i?4, 
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where the right hand sides are given by 

Ro = AiiT)X+X+W2 + A2{t){X+Wi)^ 

Ri = Bi{t)X+Wi + B2{t)WiX+Wi + B:i{T)X+W2, 

R2 = Ci(t) + C2{t)Wi + C3(t)(T^i)2 + C4t)W2 + C5{t)X+WiX-Wi, 



(4.22) 



-R3 
R4 



Bi{-t)X-Wi + B2{-t)WiX-Wi + B3{~t)X-W2, 
-Ai{-t)X^X^W2 - A2{-t){X^Wi)^, 



with T-dependent functions ^^(t), Bj{t), Ckir) which are hsted in append ix |A.3[ . These func- 
tions have been calculated from the lower order expansion coefficients (0)-(13) and from ( 4.2(]|) . 
The symmetry inherent in these expressions reflects the time-symmetry of the underlying space- 
time. 



Using the expansion ( 4. IS ) and corresponding expansions of the terms above, we decompose 

, m, satisfy the 



( 4.21 ) into the following equations. For m > 6 the coefficients (/)?_ i,, k — 



homogeneous system 



i;m,k ' 






9r4'l^,k - Wfif + ^)^lrn,k + Wf(f + 1) 



.fc= 0, 



dr(, 



3;m,fc 






0. 



"1)0? 



3;m,A; 



(1 - r)drcj>l^,, + 0|„,, - V(f + 2)(f 
The coefficients ^^.4 j., fc = 0,...,4, solve 

(1 -I- T)dr(i)l,^^k - '/'0;4,fc + 2 '/'?;4,fc = 2%/6 Ai(t) M/2;4.fe + %/6 A2(r) flfc , 



C'r</'?;4,fc + '/'?;4,fc - '/'0;4,fc + 1^/^ 4>l,i,k = 5V6S2(T)afc + VQ B:i{T) VK2;4.fc, 



9x02 



;4,fc ^ V^l;4,fc V^0;4,fc ^ 2 ^ " V^2;4,fc ~ 2 
^2;4.fc - l\^0?;4.fc + 1^ <I>IaM ^ 1^3 (t) - C5(r)] a, + C4 (t) W^2;4,fc , 
5r0i;4,fc + Cfc + Cfe - 5^0L,fc = -iV6i?2(-T) afe - V6i?3(-T) W2, 

(1 - r)9,03 4^^ + 0|4,fc - 2 0I4 fe = -2x/6 Ai(-t) W^2;4,fc - V6 A2(-r) a^, 



4,fc, 



(4.23) 



(4.24) 



with the coefficients ah defined in (4.17). The functions 0f.2 j,, A; = 0, 1, 2, satisfy 



5r0?;2,fc + 0?;2,fc + ;^</'2;2,fc = \/2Bi(t) M^i,2.fc, 

C2{T)Wxa,k. 



;2,fc ^ "'^l;2,fc ^ ^ V2-2,k 
'^■r4>2;2,k ~ yJ01;2,fe + "7J03;2,fc 



3;2,fc 



:^02;2,fe= -V2B^(-T)Wl,2,. 



^3;2,fc ^'r2;2,k 



while 02:0.0 is subject to 



5x0i;o,o = Ci{t) + [C3(t) + 2C5(t)] 6, 



(4.25) 



(4.26) 



with b as defined in ( tl.l7| ). 

These ordinary differential systems have to be integrated for the initial data ( 4.19 ) at r = 0. 
Since the equations are already quite complicated, we used t he pr ogram MapleV.4 for this purpose. 
Synthesizing the result of these integrations according to ( 4.18 ), we obtain the following concise 
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expressions for 0f on /'. 

•Pi = -(1 + t)(1 - rfX+X+Ws + ^fo{T)mX+X+W2 
+ igo(T)(X+M^i)2 + \h„{T)X+X+r^ 

<^3 ^-5(1 + r)2(l - t)^X+W3 + i/i(r) m X+W^ 

+ \gi{T)WiX+Wi + \hi{T)X+r'^ + i/ciCr) m^X+M^i, 



= -20 (1 + r)3(l - rfWa + hir) m W2 
+92{T){Wif + 3 /i2(T)r2 + k2{T) m'Wi 

+p(T)m4 + [g(T)-52(T)]6, 

= 5 (1 + r)4(l - t)2X_W3 - i/i(-T) toX_W^2 
-\gi{-T)WiX^Wi - \hi{-T)X^r'^ - \ki{-T)m'^X^Wi, 
- -(1 + t)5(1 - t)X^X^W^ + ^fo{-T)mX^X^W2 



(4.27) 



with r-dependent functions which can be found in appendix ||A.3|| . All the functions 4'f have 
polynomial dependence on r. 

The most interesting feature of this solution is its smoothness at r = ±1, which, in view of 
the singular behavior of equations ( 4.23D , (4.24) at these points, was not to be expected from the 
beginning. To explain its significance we indicate the argument which led to the asymptotic regu- 
larity condition ( 2.23| ). The Bianchi equations, which were used to obtain the evolution equations 
for the rescaled conformal Weyl spinor and, consequently, the second of the transport equations 



(4.4), form an overdetermined system. Thus there are further equations, to which we refer as to 
the constraints. In the present case the constraints take the form 



Tdr(f>l + Ux+^l - x-<j>l) - i4>l =Sl, 



Tdr4>\ 
Tdr4>l 



\{x+ 

i(x+-3 



X_03)-30: 



4-X_<^^)-3(^^ 



''S'2, 



(4.28) 



where 



Si = Fi{t)X+Wi + F2{t)WiX+Wi + Fj.{t)X+W2, 

S2^Gi{t) + G2{t)Wi + G-i{T){Wif + Gi{T)W2 + G^{t)X^WiX+Wi 

S3 = ^Fi{~t)X^Wi - F2i-T)WiX^Wi - F3{~'t)X^W2, 



(4.29) 



with functions wh ich ar e given in appendix | A.3| . As before, we obtain equations for the coefficients 
in the expansion (4.18). Together with (4.25), (4.24) these equations imply the systems 



(1 + t) (5 r2 + 3)dr(t>l,e.k + (5 r^ - 5 t2 + 5 r + 7)03. 



6,fc 



5(r- 



(1 - r) (5 r^ + 3)dr(bl,^, + (5 r^ + 5 r^ + 5 r - 7)03., - 5 (r - 






0, 
0, 



(4.30) 



and 



2(1 - rf^l^^f, + 2(1 - T)m^^, = Ti(r) a^ + T2{t) VF2;4,fc, 



4(3 + T )(1 + ■'')9r0o.4j. — ^y:L — , ) i/^0;4,fc ^ ^V-^ ^ ' ; V^4;4 

-4(3 + t2)(1 - t)9,0|4., - 2(1 + r)30|4 , + 2(1 + r)^^^^^ , = Ti(-t) a,. + T2(-t) W^2;4,fe, 

(4.31) 



with functions Ti and T2 (given in appendix | [A.3| |) derived from the functions Ri and Sj. 

It turns out that once these equations have been solved, the remaining expansion coefficients in 



(|4.18| ) can be obtained either by purely algebraic operations or by solving ODE's which are regular 
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for T e [~1, 1]. This situation is the same for all orders p > 3 in ( |4.4| ). The solutions j/(t), with 
y denoti ng in the case above the column vector with entries given by the two unknowns of ( 4.30 ) 
resp. of ( 131 ), can then be given for p > 3 in the form (suppressing here all indices) 



y{T)^X{T)X{0)-'yo + X{T) / X(r')-i 6(r') dr', 



(4.32) 



with X{t) denoting a fundamental matrix of the system of ODE's under study. The vector-valued 
function 6(t) is built from solutions which are obtained by solving the equations of lower order. 
In iQl the equations (written there in a slightly different form) have been discu sse d in general 
and the fundamental matrices X{t) have been derived. As in the case of ( 4.3(]| ), ( 4.31 ), there 
occur homogeneous as well as inhomogeneous systems for general p > 3. Thus for certain values 
of the indices (i.e. p and the indices which arise from expanding u^ in terms of the functions 
Tm * j) the functions 6(t) vanish and the solutions are of the form y(r) = X{t) X{0)~^ yo. In these 
cases some of the entries of X{t) have logarithmic singulari ties. The latter drop out of the final 
expression precisely if the asymptotic regularity conditions ( 2.23| ) are satisfied. In the remaining 
cases the entries of the matrices X{t) are polynomials in r but det(X) = c /{t) (1 — t^)p~^ with 
some constant c ^ and some polynomial /(r) satisfying |/(r)| > 1 for \t\ < 1. Furthermore, 
the column vector 5(t) has poles. However, it has no logarithmic singularities if the solutions 
of the equations of lower order have no logarithmic singularities. Assuming condition ( ^.23| ), the 
remaining potential source of singularities of u^, p > 3, at \t\ = ±1 are the integrals on the right 
hand sides of the expressions ( 4.32] ). These have not been analyzed yet. To understand the general 
situation, it is clearly of interest to study the problem for the first few values of p. Remarkably, 
in the present case, p = 3, we find that the integrand in (4.32) has poles at \t\ — ±1 and also 
outside the interval [—1,1], that the integral has poles and no logarithmic terms, but that the final 
solution is a polynomial in r. 



4.2 The detailed transformation formulae 

In this section we will determine expansions for the conformal scale factor 9 and the SL{2, C)-valued 
function A °^ w hich define the transformation from the F-gauge into the NP-gauge as described 
in section |2^|. To calculate the NP-constants in terms of the initial data we shall determine the 
values of the integrals defining these quantities by taking their limits as p — + 0. The gauge in which 
these integrals are given is based on a section C of the generators of J'^. We shall try to push this 
section to /+. The usefulness of this procedure depends, of course, on the resulting form of the 
ODE's on »/"'' which were used in |2.3| to fix the F-gauge. 

Near /+ the hypersurface J"^ can be given as the graph {r = r^ , p > 0} of the function 
T* = T^{p, t°- b) which is given by 



2n 

p 



[-DabnD''''nr- 



(4.33) 



Substituting the expansions ( [4.8| ) of Q and those of the frame vectors into the expression above, 
we get the expansion 



r" = 1 + -m p + 2 VKip2 + 0(p3). 



(4.34) 



Setting in ( 2.41 ) Z — dr , we obtain for the right hand side of this equation the expansion 



zi^Vfjevf'e) 5 



z{e) 



= -mp 



229 



24, 



-W, 



0{p% 



(4.35) 



Suppose T ~ T^ dr + T^ dp + r+ X^ + T X_ is a vector field defined near and tangent to J^ . 
Denote by T* the vector field which is induced by it on J^ . If p and i° b are used as coordinates on 
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i/"*" , one finds for T* the expression T* = T^ dp + 7+ X+ + T X- . Applying this to the gradient 
of 9 on J^ , we find that the left hand side of ( ^.4l| ) is given by 

({_2p2 + 19^p3 + 0{p'')]dp + {fX_W^ p3 + o{p^)]X+ + {f X+M^i p3 + 0(p4)}x_) (log/). 



Thus, dividing ( ^.41 ) on both sides by p^, we get a differential equation of the form r*(log/) = g 
on ^7+ with a vector field T* and a function g which extend smoothly to /+ such that T* = 
—2 dp + 0{p) near /+. For given datum /o on /+ this equation has a unique smooth solution 
which can be expanded in terms of p. As shown in our general discussion, the value of /o has to 
be constant on C to fulfill the NP-gauge conditions. We choose /o = —iryn on I^ and find for the 



solution of ( 2.41 ) the expansion 

1 



/ 



_ ■, 1 H — rn p 
2\/2 1 & ' 



191 2 ^„r 

252"^ + 5^^ 



0(p3 



(4.36) 



To obtain the matrix elements A"^ of (2.42) by using ( 2.43| ) we have to calculate the derivatives 
c*Q,(9) of the conformal factor. Using the expansion coefficients derived in pj| ], we get 



c*o,(e) = o(p4), 



c*i,(e) = V2{x+VFip3 + 0(p4)|, cV(e) = x/2{x_T4^ip3 + 0(p4)}, 
c*i,(e) = V2|-2p + 3mp2 + (8^^ - 3m2)p3 + 0{p'^)\. 



(4.37) 



Substituting these expressions into the formulae ( |2.43| ) the matrix elements A^^ and A^^ can be 
calculated expHcitly up to a U{1) phase transformation. Since the choice of the latter is not 
important for the following we choose it suitably to obtain 



A°i=p^{l-imp+(-Il^i + il|m2)p2 + o(p3)|^ A\-pf{iX+W^i+0(p)}, (4.38) 

which allows us to determine also the expansion 

El^, =V2{imp2 ^ (_^^2 ^ 21^1 )p3 + 0{p^)]dr 

+ x/2{ip2 _ 7^„^3 + (|77^2 _ ^^^)^4 _^ 0{p^)]dp (4.39) 

+V2{-|X„W^ip3 + 0(p4)|x+ + x/2{-|X+W^ip3 + 0(p4)|x_. 

To solve the differential equation for the affine parameter on the generators of J^ , we observe that 
already in the case of Minkowski space-time this parameter is a singular function of p, given by 
u° — — V2/9^^ +M°. The inspection of the expansion (4.39) suggests to search for a solution of the 
form 



V-2(-l 



7 



771 log p 



(4.40) 



This ansatz does indeed lead to a smooth regular equation for w near /+ . It allows us to calculate 
the expansion 



r 1 7 

1° = V2{ — + -777l0gp + < 

I P 3 



109 2 62^^^ , ^, 2. 



(4.41) 



where u° denotes an arbitrary constant initial datum on /+. As described in chapter |2^|, the 
matrix elements A°q and A^q can now be determined. We obtain the expansions 



A"o = pi{ix_w^i + o(p)}, x\ = p-i{- 



r^m p 



0{p')}. 



(4.42) 
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Knowing the matrix A^^ on null infinity, we can calculate the limits of the NP-spin-coefficients 
rQ^,^;^ and r^Q,QQ at /+ as p ^ 0. Substituting our expansions into the formula for the connection 
coefHcients 

^aa'bc — ^ a^ a' ^ b^ c^ff'gh ^ ^gh^ b^aa'i^ c): (4.43) 

we arrive at the expressions 

r°i,ii|,+ = lim r°^,^^ = 0, r°o,ool/+ = i™ r°o.oo = ^^- (4-44) 



The next step is to calculate the conformal scale factor 9 by solving equation (2.31). To 



determine the Ricci spinor component $22 = ^RaffEi^E^f, , we have to determine the Ricci 
tensor Raj3 of the metric g. The components of the tensor 

&a,3 ■■= ^R(a0) - -^QapR + ^-R[a/3] (4-45) 

in the frame {c*^/}, where Rap resp. R denote the Ricci tensor and the curvature scalar induced by 
the Weyl connection V with coefficients T^^ — T^*^^ + (5f /^ + 5^ fa — ga-tf^ (cf. [g|) , are among the 
variables of the conformal field equations. Thus they are known to 3''''-order in the p-coordinate. 
From the general transformation law 

Rap = Rap - 2V(,//3) + 2fjp - gapi-^-^P + 2PP) + 4V[„/^], (4.46) 

we get the relation 

Qap = 2 (^"/3 - 6^"'^^) ~ ^'^^" ^ ■^"•^'^ ~ l^9o,pf-tP- (4.47) 

From this we derive the expression 

<i>22 = QapEl^El^, + Ely {El^P) - (El^Pf. (4.48) 



Substituting here (|4.39| ) and the expansion of the one-form / obtained from the solution of the 



field equations we get the expansion 

*22 = \mp^ + (-^^' + y W^i) P' + 0{p^) (4.49) 

on J^+. 

On /+ is induced in our gauge the standard S'^-metric. Therefore we solve equation ( 2.31 ) with 
the initial condition 

lim 6 = 1. (4.50) 

For the conformal scale factor we obtain then the expansion 

9 = 1 + \mp + {^W, + i|m2) p2 + o{p% (4.51) 

By the choice of the initial value for the conformal factor the scale function p appearing in the 
gauge transformations is also fixed with 

p=l on J+. (4.52) 

In the conformal gauge characterized by the conformal factor 6* := 9Q the generators of null 
infinity are expansion free. Proceeding as indicated before, we construct the NP-frame {-E'^/}. 
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Observing the expansions ( p.42|) and ( |2.45[) of the null vectors E^y resp. E'^^, and taking into 
account the properties of the conformal rescahng we get the relations 



A°i 



3-i \0 „ic 



A"ie'S A\^e-2X\e 



(4.53) 



with function c, characterizing the phase freedom, which will be fixed later. Using ( 4.3g ) and ( 4.51 ) 
we get the expansions 

A\ = pi {l - Imp + {§ m^ ^ 2Wi) p^ + 0{p^)^e"', A\ = pi|i X+W^i + 0(p)}e*^ (4.54) 

from which we derive in turn the expansion 

Ely ^V2[\mp^ + (-m2 + 2Wi)p'^ + ©(p^)}^, 

+V2{y^ - 2mp3 + (^m2 - fWi)p^ + Oip^)]dp (4.55) 

+V2{-^X_Wip^ + Oip^)}x+ + V2{~IX+Wip^ + 0{p'')}x_, 

of the vector field E^y tangent to the null generators oi J"^ . Furthermore the new affine parameter 
has the form 



= V2|-i-™-^'-' ^'''-^' 74,,, ,„_,^, 



4mlogp + <+ ( —m' + —Wi]p + Oip')}, 



with a free constant u*. Using the formula analogous to (2.44) we derive 



A\=pi{-^X_W,p + 0{p^)y 



Ai 



'{- 



jmp 



0(p^)}c 



(4.56) 



(4.57) 



To determine of the phase factor e^*'^ we solve equation ( 2.33 ) along the generators of null infinity. 
Expanding the right hand side, we get 

Ely{c)^ 23m {A\l^\,A\A\r}f,^„~A\Ely{A\)+A\El,,{A\)}, (4.58) 

where A°j, has been obtained from the matrix A°j above by setting c — 0. Substituting the 
known data into the equation above, the solution c which is needed to satisfy the gauge condition 
-'^11'oiIj ~ 0' i^ found to have an expansion 



which entails the expansions 

e- = l + 0(p2), i?*^,(e-) = 0(p3), i?oV(e'^) = 0(p2). 



(4.59) 



(4.60) 



The matrix elements A"j, are now d eter mined on null infinity to the precision needed in our 
later calculations, but in the definition ( |3.5| ) of the NP-constants appear some o f the transversal 
derivatives -Eqq/(A°j) of the matrix elements as well. Using the general formulae ( 2.47 ) we get the 
expansions 



'oV(A"o) - V2pi{'-§X^W, + 0{p)}, E',,,iA\) = V2p-^{l + §mp + Oip^)}, 

oV(A"i) ^ V2p--^{\ + %mp + 0{p^)}, i?oV(A\) = V2p'^{-%X+W,+0{p)}, 



(4.61) 



where we have taken the expressions ( fl.60| ) for the phase factor into account. 

The transversal derivative of the conformal scale factor E'^, (9) is fixed on null infinity by the 
requirement i?[g'^]|j+ — 0. Thus it has to satisfy equation (2.3S) with initial datum 



E'„,,{9)\j^^lunep-'r 



1-po 

lo'oo ~ "^^A^ lo'oo- 



lim r? 



(4.62) 
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Given the matrix A*^^ and the conformal scale factor 9, all the terms appearing in equation ( |2.38[) 



can be calculated in a straightforward way, with the exceptio n of t he curvature scalar R[g], whose 



calculation requires some explanation. Contracting equation ( [4.47| ) we get the identity 

R[g] - 6 iOaa'bb' + ^aa'fw + faa' fw)e^''e^'''' , (4.63) 

where 

^aa'fbb' = Caa'ifbb') - {^la' cb^b' c' + ^*aa'c'b'<^bc) f"" ■ 

Expanding these quantities we get 

/ ^ X (4.64) 

F* = [§ m2 - 4Wi + la-X+Wi - ^a+X.Wijp^ + 0{p^), 



which entail with ( ^.38 ) the expansion 



^0*0' {e) = V2l—m+(—m^-4Wi + -a^X+Wi - -a+X^wA p + 0{p^) | . (4.65) 

Given the expansion above, we can calculate expansions of various quantities of physical interest, 
such as the Bondi energy momentum, the angular momentum, and the radiation field on J'"'". Since 
the coefficients in these expansions are given directly in terms of the initial data on the Cauchy 
hypersurface S, the expansions contain information about the evolution of the field over an infinite 
range. As an example we will calculate below the NP-constants. 

We close this section with a remark on the BMS group, the group of transformation between 

different Bondi-type systems. It was shown in |Q that for solutions for which the the condition 

lim rr 1 Q,,Qg — could be realized at space-like infinity, where the subscript "e" is to denote 

the electric part of the considered spin-coefficient, one can single out the inhomogeneous Lorentz 
group as the group of transformations preserving this condition. It turns out that under our 
assumptions, which include in particular the time-symmetry of the solution, the even stronger 
condition lim FJ^^/qq = is satisfied. This means that for our solutions there is a natural way 

u* — > — oo 

to single out the inhomogeneous Lorentz group as asymptotic symmetry group. 

4.3 The NP-constants in time symmetric space-times 

Using the formulae of the previous chapters we can express the NP-constants in terms of the initial 
data for the corresponding time symmetric solutions. All the quantities appearing in the integral 
(|3.5|) are known in terms of the initial data to the precision needed to perform the limit p ^ 0. 

W e have to express the spin-2 spherical harmonics 2^2,^1 in terms of the functions T^ ^. By 
( |3.3|) the defi nitio n of the 9-operator is based on the choice of the complex null vector field E^y . 



In appendix [1A.1| we have appHed the standard choice and derived the relations between the 
operators X+ and 8 and between the spin-2 spherical harmonics 2^2, m and the functions T^ j.. 
By this choice we should have E'y — -y^X^ on /+. However, calculating the vector E'^ in the 
conventions used above, we get 

SoV|/+ = ^X_. (4.66) 

There are two causes of the difference. We fixed the phase factor such as to simplify the calculations 
and the conventions used in the F-gauge and the NP-gauge are such tha t one has to swap the two 



spinors of the dyad to get from one to the other convention. The form (|4.66D o f E'y corresponds 



to — i-\/2m, if m denotes the standard complex null vector used in appendix \AA\. This means that 
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( 4.66| ) co rresp onds to the operator — i8 instead of 5 discussed in the appendix. Observing this and 



(|A.g| ) in (|3.4|) we obtain the formula 

G™ = j2-m(5^)i J Tf-f i?oV(<^o)M for m = -2, . . . , 2, (4.67) 



where /i = j^dtS da is the Haar-measure on SU{2). 

To calculate ( 4.67| ) we expand the integrand in terms of p and take the limit as p ^ 0. For 



this we have to determine for i?QQ/(0o) only the terms of order 0(1). In the limit only these terms 
give a contribution while the terms of order p^^ cancel each other. Using the explicit results of the 
previous chapters we arrive after some lengthy but straightforward calculations at the expression 

/-o „/ 5 o 635 1905, ,, 16 \ (4-68) 

f^tT[-^X-X-r^ + -^mX_X_W2 - ^- (X.Wi)^ + —X^X^W^jp. 

Expanding the functions in the brackets in terms of the functions T^, and using the orthogonality 
relations satisfied by these functions we can p erfo rm the integr ation . All terms expect the last one 
give some contributions. Using the formulae (4/7), ( 4.12 ) and ( 4.16|) we get the final expression 



Gm\j+= ^-(157r)5 j 127 (m VK2;4.2-rn " 6 fls-m) " "^ i?2-m [ > (4-69) 

where the coefficients a2-m, which are quadratic in Wi;2,fc, are given by ( [4. 17] ). We note that the 
structure of this more general expression is essentially the same as that of the expression obtained 
by Newman and Penrose in the case of static and stationary solutions. 

5 Concluding remarks 

We have seen that, under the assumptions explained above, certain fields which are given near 
space-Hke infinity in terms of Bondi-type systems can be expressed in a straightforward way in 
terms of the gauge conditions used in ^ and can thus be related directly to the structure of the 
Cauchy data which give rise to the space-times by Einstein evolution. The calculations involved 
are quite lengthy but taking into account that we relate quantities which are obtained by a non- 
linear evolution over an infinite domain of space-time to the data from which they arise, the overall 
structure of the argument is surprisingly simple. 
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A Appendix 

A.l X+ and the 5-operator 

In this section we describe the relation between the operators 9, 5 introduced in [|l3| and the 
operators X+, X_, X used in |Q]. 

Consider on the group SU{2), which is diffeomorphic to S^, coordinates {x, y, a} such that 
outside a set of measure zero the general group element t°^ £ SU{2) is given by 



1 / e*"_ ie~'°'C 



^b^ r. TT 1 ^r^ioL/- „-ia ] J (A.l) 
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with C, = x + iy. Then a is a parameter and x and y are constant on the orbits of the the subgroup 
C/(l). The tangent vectors dx, dy, respectively da at the unit element coincide with the generators 
ui, U2, and U3 of the Lie algebra of SU{2). Writing P = 5(1 + CO) we get for the corresponding 
left invariant vector fields the expressions 

Zui = P cos{2a)dx + P sin{2a)dy + ^[xsin{2a) — 2/cos(2a)]9Q, 

Zy,^ = -Psm{2a)dx + P cos{2a)dy + ^[ysm{2a) + xcos{2a)]da, (A. 2) 



7 — Lf) 



whence 



X+ = -Zy, - iZu, = e2-{-jV2(TO - ^C^a}, X = -2^Z„3 = -id^, 
X- = -Zu, + iZu, = e-2»{jV2(m + ^Cd^}, 



(A.3) 



where the vectors m = a/2F9^ and m = \/2P9^ define a complex dyad tangent to the surfaces 
{a = const.} which is null with respect to the standard S'^-metric ds"^ = P^^d(d( on these 
surfaces. 

We may identify SU{2) with the spin frame bundle over the base manifold S"^ with struc- 
ture group U{\). The section {a — 0} can be identified with the base manifold (with a point 
omitted). Here we take the complex null frame {m,Tfi} defined above, where a group element 
u'\ = diagie^" ,e~'^") G U{1) acts as u{{m,fh}) = {e^'^m, e~^'"m}. A function ry on S^ is said to 
have spin weight N, if it can be decomposed as r]\^^a = e^^*"?7o, where the function 770 is indepen- 
dent of the parameter a along the fibers. The 9-operator is defined by the complex null vector m 
and acts on a spin-N function as 

Qrj\ca = V2{m{rj„) + Nrj„ m''m''6pm^}e^(''+^'>"' - V2{m{7jo) + ^iVC^o}e'^^+'^", (A.4) 

v2 

where S denotes the Levi-Civita differential operator induced by the standard 5'^-metric. This 
means that dr] has spin weight TV -I- 1. (This treatment of the functions with spin weight and the 5 
operator is a bit different from the one which can be found in the literature (cf. ||l^, g, 0), where 
the expressions are evaluated on some cross-section of 5^.) 

The horizontal lift of the vector m defined with respect to the Levi-Civita connection S is given 

by 

"Ifflca =™" ^C9a. (A. 5) 

This means that the 8-operator on S^ is given by 

9|^ „ = V2e2^"mff. (A.6) 



Comparing the formulae (A.3), (A. 5) and (A^) we get the relations 



X+ = -iQ, X^^iQ, a: = -[9,9]. (A. 7) 

The spherical harmonics Yi^m are defined as an orthogonal function system on the sphere S'^. 
They can be extended to S^ as functions with zero spin weight, i.e. they became independent on 
the parameter along the fibers. This means that they can be expanded as Y/.m = /.^CkjT^j/ /^ in 

terms of the functions PJ/.- The spherical harmonics satisfy the equation dQYi^rn — —l{l + l)^,m, 



so using the relations (A. 7) and (2.22) we arrive at the relation 
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Taking into account the explicit coordinate expressions of the group elements one could determine 
the ex pa nsio n coeffi cient s Cj. Using the definition of the spin harmonics sYi^m (cf. ^) and equations 



( 2.22 ), (A.?) and (A^) one can also derive the relation between the functions sYi^m and the 
functions T^^ j,. We shall only need the transformation formulae 

\^ / ■\4:—m,t 5 \2rp2—m 

^2,m- (-«) (i^j J4 2 : 

V — ( 4\2-'m( 5 \2rp2-m -i^ _ / ■\2-m ( 5 \2rp2-m 

2J'2,m- (-«j i47J J4 : -2J'2,m - (-^j (,4;;pj J4 4 ■ 

A. 2 Some useful spinor identities 

Here we describe irreducible decompositions of spinors with four unprimed indices in terms of the 
"primary spinors" e^bcd' ^abcd, Xab, Vab, Zab and eab, where 

r , - A/2f °f ^ VI,- Lf If 1 z ^ - ^f °f ° 

•i'ab — V ^t(-^ tj^^ , L/ab — ^'^a '^b ' ^ab — ^'^a "^6 ' 

i _ {e f g h)i , _ \ ■ J 

^abcd — ^{a ^b ^c ^d) ' "ohcd — —<^a(c^d)b- 



It is well known that a spinor Aabcd satisfying Aabcd — A(^nb)[cd) = —^cdab can be decomposed in 
the form Aabcd = ^acAbd + ^bdAac with Aab = 5 ^a/b''^ = ^(ab) ^^d that a spinor Sabcd satisfying 
Sabcd = S'(ah)(cd) = Scdab Can bc Written in the form Sabcd = Si^abcd) + \habcdS with S := S'^^''-^. 
It follows from this that an arbitrary four index spinor with symmetries Xabcd — X^ab){cd) can be 
expanded in terms of ej,^^^, eacXbd + f-bdXac, f^acUbd + f-bdVac, eacZbd + f-bdZac and habcd- 
The following relations were frequently used in the calculations: 

UabXcd = ~^abcd ~ r-K^acVbd + ^bdVac), Zab^cd = ^ abed + r-\(-acZbd + ^bdZac)] 

XabX'''' = -l, Xahy"'' = 0, a:a6Z°' = 0, yaby"' = 0, y,6z'^'' = -i, Z.fc^'^'' - 0; 
^a ^b/ = 2^''b7 Va ^bf = -^Uabi Z^ Xbf = —-^Zab: 

c-O ^cd _ n J.0 ,,cd _r , e-0 „cd n ^1 cd — _ 1 ^ , 

'^abcdX — "i ^abcdy — "^ab, '^abcd^ " ^' ^ abcd"^ — 2'^«b, 

^abcdy ^ ^Xab, ^ abcd'^ — "' ^abcd-^ — z^ab, fcafccdS/ " glJab, 

J.2 cd — 1^ , £.3 r^cd 1,1 , e-3 ,,cd n ^.3 „cd _ 1^ , 

^abcd-^ — ^Zab, ^abcd-^ — 2^°-^, ^ abcdti " ^' ^ abcd'^ — A^ab, 

^abcd^ — 0, ^abcdV — 0' ^ohcd^ = ^Vab] 

X(abXcd) = ^^abcd^ X(^abycd) = ~^abcd' X(^abZcd) = ^ abcd^ 

14 12 10. 

ViabVcd) ~ 2^abcd' V(abZcd) — ~2^abcd' Z^abZcd) — 2^abcd' 

^{a ^b)cdf~^^abcd^ ^ (a ^ b)cdf ~ ^^^"-^Xcd^ X ^^^ ^b)cdf = ui^acXbd + ^bdXac) , 

^{a ^b)cd/ ~ 275^"''^'='^' ^{a ^b}cdf ~ ~ ^^ abcd^ '''ab{c ^d)f=2\^acXbd+^bdXac)] 



y^d ^c)abf ~ 2V2 ^abcd + 24 V^acJ/bd + Cbdyacj J Z^^ ^c)abf~ 2V2^abcd "I" 24 V^"c^bd + Cbd^acj ! 
^ ab ^cdef~~12^abcd^^^y^acZbd + ^bdZac), ^ ab ^cde/ " ~ 12 ^abcd + ^^f '^"'=^'"' "*" ^^'^^"'^^ ' 



^2 ^2 obcd _ j_ J.2 e/ 2 __lc2 J- J_h 

^abcd^ ~ 6' ^ab ^cdef~ 6 ^ abed 1" 18 "bed- 
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A. 3 The detailed expressions for u^, p = 0, 

The T-dependent functions occuring in (|4.2|). 



c'^^(t) —■m{ 



4 ^3 



3 ' >' 



.±1 



(r)=m(r2-ir4), ^^(r) 



^r^), 



ii:i(T)=TO(-12T + 4T3), Fi(T) = imT4, t^ir) 

ri(r)=6™(l-r2), </.i(r) = -12(l-r)2, 0i(r) 

(?!)^(r) =-36 + 36r2. 

The r-dependent functions occuring in (|4.3|). 



= \/2TO(iT 4 

= \/24tto, 

= -to2 (18^2 _ 3^4 



.±2 



„±2 



5| 



Fi 



T? 



n 



m2 (-2x3-3x5 + 1 r7 - 


1 , 
7 


r^), 


8r3-Zr5-|r^ 






rn2(-2T2 + 3r4-|r*^ + 


1 

14 


r«) 


-6r2-ir4 + A,6^ 






V2(6r2-|r4 + ^^6)^ 






-36T2 + llr'' + |T^ 






144t + 72t3 ifrS, 






-V22r3, 






m2(-2T2 + ir4-|T6 + 


1 
7 


t'\ 


3r4 + |T6, 






y2(48T-16T3), 







= n? (-12 + 12 t2 - IOt^ + I T^), 
= m2(4r2-|T4 + |T6), 
= 24-12r4, 
= 4m(^ 



'37 T-6 

-10 ' 



41 5 
5 



41 4 
2 ' 



46r3-18T2), 



= 6(-Ar8 + Mr6-15T4 + 6T2)m3, 
= -72(l + T)2(-l + r)2. 



^02 



(r 



ci2(r 

C2^^(r 
5?(r 

Sl{r 

mr\ 

KKr 

FKr 

tl{r 

tl{r 
Tiir 
Tl{r 



16r3-f rS + fr^, 



= TO(-4T2 + |r4), 



.4 _L 3 ^6 
5 



V2m^(|r4-|r6-^r«), 
V2(-|r4 + 3r2-|,r6), 
TO2(24T-8r3 + 4r5- ^ r^) 



m^-fr 



20 ^3 I 8 ^5 20 ^7 



20 7) 



= 2t 



;r+f r5, 

4 6 _6 

5 ' 



x/2m2(-12T-fT3 + |T5), 
\/2(24r + 8r3), 
72-72t2 + 36t4, 
-\/26t2, 

-(-l+r)^ 

16(1 + t)(-1+t)3, 
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The r-dependent functions occuring in (|4.2(]| ) 



cf 


t) 


cf 


t) 


„03 
C3 


t) 


„03 

C4 


t) 


cf 


t) 


cf 


t) 


cf 


t) 


43 


t) 


cf 


t) 


cf 


t) 


„±3 
C3 


t) 


cf 


t) 


cf 


t) 


cf 


t) 


5? 


t) 


5l 


t) 


5l 


t) 


5I 


t) 


SI 


t) 


SI 


t) 


5| 


t) 


^1 


t) 


K? 


t) 


Kl 


t) 


Kl 


t) 


Kl 


t) 


Kl 


t) 


Kl 


t) 


Kl 


t) 


Kl 


t) 


Kl 


t) 


Kfo 


t) 



= (3t3 + 18t5 
= (-44 r3 



283 7 _ 1510 9 I 2972 11 
21 189 ~'~ 2079 



693 



^)'m? 



5|8^5^2M^7_^^9 + |^11)^^ 



;r3-2§r5 + 4^T7~^r9, 



(-20r3-6r5 + ^r 



-3 _ 6 -5 
= (12t2 + 15t4 



439 ^7 _ 573 9 

280 ' 



ij^")m, 



4 ^7 ^ 4 9 

7 ' ^ 7 ' ' 



+ ^T^)vn?, 



= -72T2 + 18T4-f t6 
= -36r2 + 3T4^^-6 

= (18r2 + 12T4 
= -12r2 



31 ^6 I 3 8 
5 ' '^ 2 ' 



■T^°)m, 






— r^ T-2 _ 33 T-4 1 50 6 „ 515 ^8 I 25 _10 _ J3£ _12^ _,3 

— {2 ' 2 ' ""^ 3 84 ' "'' 27 693 '' ' 



(-4»r 

36t2-12t4 

-3r2-2r4 



2^105r4-^T6 + ^r«-|TiO)m, 



24 _6 _ 6 _8 
5 ' 7 ' ' 



3 „6 
5 



14 ' 



= -9T2-2r4 + f T^ + ^L^s 



(108 t2- 168T'^ + 86r' 



6 _ 39 _8 
5 ' 



■ T^'^) m, 



= -72t2 



(- 



72 T-6 _ 4 8 
5 ^ 7 ^ ' 



9 ^2 37 ^4 I 19 ^6 827 , 



168 



355 10 
378 



^ri2)y2 7 



77 



69 4 _ 333 6 1 1999 -8 , 13 10^ 
2 ' 20 ' 560 "•" 80 ^ 



V2 



= (6t2 

= (18T2-6T4 + f T6-fT8)x/2, 



TO, 



(6t2-2t4 ' 2.8 



r")\/2m, 



-6r 



S^3_^M^5 + 4^7^ 



(144t+12t3 



351 .5 I 237 .7 17 .9 



T ) TO, 



= -96 r + 16 r3 + f t5 - M T-7 ^ 



(-54r + 12r3-216r5 



796 .7 440 .9 1 16 .11 



21 



+ ifr")TO3, 



864 5 I 288 .7 



= (576t - 216r3 + if2 r^ - 114 ^7 + f r9) TO, 
= -432T + 288r3. 

= (40 r3 - 16 t5 + ^ t7 - i|^ r'J + J^ r" ) TO^ 
= (-240 r^ 



582 .5 _ 218 .7 _L 23 .9 
5 



r^-^T' + fT^)m, 



= (-4t3 + |t5)V2, 
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Fiir 
Fiir 
Fiir 
Fiir 

t\[T 

tUr 

tUr 

tlir 

tUr 

THr 

TUr 

Tiir 

nir 

Tiir 
Tiir 
THr 
Tiir 
Tiir 
Tioir 



9t2 + 2t4 



3 ' 



26 _8 
7 ' 



20 _10 

21 ' 



^T^^)m^ 



693 



-60r2 + 36T4-12r6 + ^r8-fri0)m, 



24 6 
5 ' 



12 8 

7 ' ' 



-12t'-6t^ 



7 6 , 169 8 , J_ 10^ 

2 ' ^ 56 ' ^ 40 ' '' '"■' 



144 ^5 
5 



r-)V2, 



294 ,5 
5 ' 



§r^)V2r 



36t + 20t3 + 46r^ - ^ T^ + fg r^) V2m^, 
-312 r - 24 t3 - f r^ - f r^) \/2 m, 
144T-96r3- 
-96 r- 12 t3 
48T-f t5)V2, 
3 + 9r2-3r4-r6, 
-72 - 36 r2 + 81 t4 - 423 ^6 + ^ ^8) ^^ 

=8-24T2 + 16r^ 

27- 18r2 + 180r4- 134 t^ 



204 8 
7 ' 



f ri")m3, 



-288 + 216 t2 - 558 r4 + lf6 r6 - Ml ^8) m, 
216 - 216 r2 + 216 t" - 72 r^ 
-24r2 + 16T4-f r6 + If r8- 4 Ti0)m3, 

144 t2 - 102 r4 + i|8 ^6 _ |7 ^8) ^^ 



z / T _j r 20 ^ 
-12r2 + 6T4)\/2. 



|.r«)V2m, 



The T-dependent functions occuring in ( |4.22| ). 



Ai(t 
A2ir 
Bi{t 

B2ir 
Bsir 

Ciir 
C2{r 
C^ir 
Ciir 
C,ir 



97 _4 
2 ' 



6 T-5 
5 



169 ^6 
5 

2592 ^5 



= (36 r - 78 r2 + 82 r^ 

= -648 T + 1728 t2 - 1692 r^ + 432 t^ 

== (108 r - 234 r^ - 396 r^ + 1503 t* - 579 t^ - 

I 11682 7 I 40413 ^8 _ 2591 _9 , 177 _-10\™2 
I Qc; ' ~r r.p.rt ' 70 SO -^ ' 



2M^r_^^^8)^^ 



2286 6 
5 



162 



Z56^7_ . 
5 5 



20 



35 



560 



-648 T -I- 1404 r^ - 540 r^ - 810 r"* -f if^ r^ -I- if^ r^ - 108 r^ 



= (-72 T -f- 168 r2 -f 24 r3 - 274 r" -f 120 r^ -f 226 y 
= (-27r-K342T3-696T5 



32 r^ 



2598 _7 _ 4555 ^9 
7 63 



1079 ^11 
231 



) to"* 



= -1296 T + 2376 t3- 4^r5-K216T^ 
= (-432 T -h 792 T^ 



3072 ^5 
5 



^r7)m, 



= -216t + 1Q8t^ + ^^ t"" + ^ t' . 



324 ^7 



108 _8 
5 ' ' 



I m, 



REFERENCES 34 



The r-dependent functions occuring in (|4.27| ) 



/o(t) = -18 + 216 T^ - 240 t3 + 18 r"* - 48 t^ + 204 r^ - 144 t'^ + 30 r^, 
/i(t) - -9 - 216 r2 + 696 r^ - 198 r^ - ^Mi ^5 + 9|4 ^6 ^ 9|6 ^7 _ 4ii ^8^ 

/2(t) =-3 - 216t2 + 372r4 - 9f6 t6 + 219 rS, 

5o(t) - 108 - 1944 t2 + 4752 r^ - 5724 ^^ + m r^ - ^ r^ + ^M ^7 _ i|8 ^8^ 

51 (r) = 54 - 972 t^ + 1620 t^ + 378 r* - iif^s r^ + ^ r^ + if r^ - if r^, 

52(t) = 18 - 540 r2 + 972 r^ - 2|8 r^ + if r^, 

/io(t)= i, /ii(t)= |, h2{r)^ i, 

fcl(T) = 108 t2 - 276r3 - 129t4 + ifZ r5 - 3289 ^6 _ 9439 ^7 + 32803 ^8 + ^ ^9 _ ^ ^10 

fc2(T) -252t2 _ 942t4 + 3^t6 - 6||ir8 + f r 
p{t) =-f r2 + iIiT4-116r6 + i||2 

g(T) =216^8_ 576^6_^g48^4„8g4^ 



9439 
35 


T^ + 


32803 
280 


r« 


1079 
2772 


rl^ 







The T-dependent functions occuring in ( 4.29| ) 



1071 _4 I 4077,5 2639,6 18878,7 i 113287,8 i 1389 9 15087 IQ-) ^2 



F-. (t') - r72 -7-2 - 1211 T-4 I 4077 5 2639 ,6 18878 ,7 , 113287 ,; 
riyT) — ((/.T 2^+5^ 20 ^ 35 ^ "I" 560 ^ 

F2(t) = -864 t2 + 1584 t3 - 810 t^ - if^ r^ + 8|2 ^6 + 4|2 ^7 _ i^s ^8 
F3(r) = (-36T2-40r3 + 156T4- 8f T^ + lf r6 + 4|6^7_ 198^8 
Gi(r) = (f r2 + i|i r* - 348 r6 + 6||5 t8 - i^ r^" + ff ri2) m\ 
G2 (r) = (-144 t2 - 1071 T* + 3^ r6 - 2||p r^ + 24|9 ^lO) ^2^ 
G3(r) =1116x4-468x6+648^8, 
Giir) =174mT4-lf4mr6+6fimT8, 
G5(r) =432x2-234x4+306^6+216^8. 



20 ' 560 



m. 



The r-dependent functions occuring in (4.31) 



Ti(t) = -46656 r+ 124416 r2- 138240x3 + 51840x4 + 31104x5- l^^pe ^6 



+ 12096T^-^T8 + i2|68r9^ 

r2(T) = (5184 T - 3456 t^ - 5088 m r3 - 6048 r^ + 12288 r^ + 384 r^ - 4128 r^ 
+ 1824 r8- 1344r9 + 384TiO)m. 

References 

[1] Bondi, H.,van der Burg, M. G. J. and Metzner, A. W. K. 1962 gravitational waves in general 
relativity, VII. Waves from isolated systems.. Proc. Roy. Soc. Lond. A 269, 21-52. 

[2] Dain, S. In preparation. 

[3] Friedrich, H. 1995 Einstein equations and conformal structure: existence of anti-de Sitter-type 
space-times. J. Geom. Phys. 17, 125-184. 

[4] Friedrich, H. 1998 Gravitational fields near space-like and null infinity. J. Geom. Phys. 24, 
83-163. 



REFERENCES 35 



[5] Friedrich, H. 1998 Einstein equations and conformal structure. In: The geometric universe: 
science, geometry and the work of Roger Penrose, S. A. Hugget et al. (eds.), Oxford University 
Press, Oxford. 

[6] Friedrich, H. 1998 Einstein's equation and geometric asymptotics. In: Gravitation and Rela- 
tivity: At the turn of the Millennium, Proceedings of the GR-15 Conference, N. Dadhich, J. 
Narlikar (eds.), lUCAA, Pune, India. 

[7] Friedrich, H. and Schmidt, B. G. 1987 Conformal geodesies in general relativity. Proc. Roy. 
Soc. Lond. A 414, 171-195. 

[8] Goldberg, J. N. 1967 Invariant transformations and Newman-Penrose constants. J. Math. 
Phys. 8, 2161-2166. 

[9] Goldberg, J. N. et. al. 1967 Spin-s spherical harmonics and 8. J. Math. Phys. 8, 2155-2161. 

[10] Ko, M., Newman, E. T. and Tod, K. P. 1977 H-space and null infinity. In: Asymptotic structure 
of space-time, F. P. Esposito and L. Witten (eds.). Plenum Press, New York. 

[11] Newman, E. & Penrose, R. 1962 An approach to gravitational radiation by a method of spin 
coefficients. J. Math. Phys. 3, 566-578 & 4, 988. 

[12] Newman, E. T. and Penrose, R. 1965 10 exact gravitationally-conserved quantities. Phys. Rev. 
Lett. 15, 231-233. 

[13] Newman, E. T. and Penrose, R. 1966 Note on the Bondi-Metzner-Sachs group. J. Math. Phys. 
7, 863-870. 

[14] Newman, E. T. and Penrose, R. 1968 New conservation laws for zero rest-mass fields in 
asymptotically flat space-time. Proc. Roy. Soc. Lond. A 305, 175-204. 

[15] Newman, E. T. and Unti, T. W. J. 1962 Behavior of asymptotically flat empty spaces. J. 
Math. Phys. 3, 891-901. 

[16] Penrose, R. 1963 Asymptotic properties of fields and space-time. Phys. Rev. Lett., 10, 66-68. 

[17] Penrose, R. 1965 Zero rest-mass flelds including gravitation: asymptotic behavior.. Proc. Roy. 
Soc. Lond. A 284, 159-203. 

[18] Penrose, R. 1967 Conserved quantities and conformal structure in general relativity. In: Rel- 
ativity and astrophysics, J. Ehlers (ed.). Am. Math. Society. 

[19] Penrose, R. 1982 Quasi-local mass and angular momentum in general relativity. Proc. Roy. 
Soc. Lond. A381, 53-63. 

[20] Press, W. H. and Bardeen, J. M. 1971 Non-conservation of the Newman-Penrose conserved 
quantities. Phys. Rev. Lett. 27, 1303-1306. 

[21] Robinson, D. C. 1969 Conserved quantities of Newman and Penrose. J. Math. Phys. 10, 
1745-1753. 

[22] Sachs, R. K. 1962 Gravitational waves in general relativity, VIII. Waves in asymptotically flat 
space-time. Proc. Roy. Soc. Lond. A 270, 103-126. 



